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This  research  presents  a  new  mathematical  approach  for  finding  an  optimum 

solution  to  a  singular  inverse  radiation  transport  problem  involving  the  nondestructive 

assay  of  radionuclide  inventory  in  a  nuclear  waste  drum.  The  method  introduced  is  an 

optimization  scheme  based  on  performing  a  large  number  of  numerical  simulations  which 

account  for  the  counting  statistics,  the  nonuniformity  of  source  distribution,  and  the 

heterogeneous  density  of  the  self-absorbing  medium  inside  the  waste  drum.  The 

simulation  model  uses  forward  projection  and  backward  reconstruction  algorithms.  The 

forward  projection  algorithm  uses  randomly  selected  source  distribution  and  a  first  flight 

kernel  method  to  calculate  external  detector  responses.  The  backward  reconstruction 

algorithm  uses  the  conjugate  gradient  with  nonnegative  constraint  or  maximum 

likelihood  expectation  maximum  method  to  reconstruct  the  source  distribution  based  on 


xn 


calculated  detector  responses.  Total  source  activity  is  determined  by  summing  the 
reconstructed  activity  of  each  computational  grid.  By  conducting  10,000  numerical 
simulations,  the  error  bound  and  the  associated  confidence  level  for  the  prediction  of  total 
source  activity  is  determined.  The  superiority  of  using  multiple  energy  peaks  in 
backward  reconstruction  algorithms  is  analyzed. 

The  accuracy  and  reliability  of  the  simulation  model  are  verified  by  performing  a 
series  of  experiments  in  a  208  t  waste  barrel.  A  number  of  Ir''^  seeds  are  used  to  simulate 
unknown  distributions  of  the  source  in  the  barrel.  Density  heterogeneity  is  simulated  by 
using  different  materials  distributed  in  37  egg-crate  type  compartments  simulating  a 
vertical  segment  of  the  barrel.  Four  orthogonal  detector  positions  are  used  to  measure  the 
emerging  radiation  field  from  the  distributed  source.  Results  of  the  experiments  are  in 
full  agreement  with  the  estimated  error  and  the  confidence  level  which  are  predicted  by 
the  simulation  model. 

This  research  presents  a  way  of  defining  the  error  bound  and  its  associated 
confidence  level,  which  is  inadvertently  ignored  by  the  traditional  waste  drum 
characterization  method.  The  concept  of  using  multiple  energy  peaks  with  advanced 
backward  reconstruction  algorithms  also  provides  a  more  accurate  prediction  ability  than 
does  the  traditional  method. 
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CHAPTER  1 
INTRODUCTION 

The  inverse  radiation  transport  problem  discussed  in  this  dissertation  is  to 
characterize  radiation  sources  inside  a  black  box  by  using  remote  sensing  methods.  This 
problem  exists  in  such  fields  as  in  vivo  radionuclide  assays,  remote  sensing  for  geology 
or  astronomy,  and  chaos  studies  for  heat  conduction.  In  nuclear  engineering,  the  most 
prominent  task  related  to  the  inverse  radiation  transport  problem  is  the  characterization  of 
radioactive  sources  inside  a  nuclear  waste  barrel.  Before  shipping  the  nuclear  waste 
barrel  to  disposal  sites,  barrel  contents  need  to  be  determined  as  high  or  low  level  waste 
based  on  the  types  and  the  amount  of  radionuclides  inside.  The  cost  of  processing  high 
level  waste  is  usually  much  greater  than  that  of  low  level  waste.  Therefore,  accurate 
characterization  of  the  waste  barrels  is  not  only  essential  environmentally  but  also 
important  economically  for  the  nuclear  industry. 

Most  nuclear  waste,  including  commercial  power  plant  waste  and  defense  waste, 
is  stored  in  cylindrical  steel  drums  (208  I  waste  barrel)  about  57  cm  in  diameter  and  88 
cm  m  height.'  Millions  of  waste  barrels  are  waiting  to  be  characterized  in  the  United 
States.  Therefore,  gamma  ray  detection  time,  speed  of  analysis,  and  accuracy  of 
evaluation  are  very  important  factors  in  satisfying  the  regulatory  requirements  in  a  timely 
manner.  For  most  applications,  using  far  field  measurement  and  assuming  the  drum  as  a 
hne  or  a  point  source  are  not  practical  because  of  the  long  measurement  time  associated 


with  the  long  measurement  distance.  For  a  near  field  measurement,  the  estimation  error 
for  the  total  source  activity  is  usually  very  large  because  the  detection  rate  is  very 
sensitive  to  the  source  distribution  and  the  density  of  the  self-absorbing  medium.  Figure 
1-1  shows  a  generic  configuration  of  a  distributed  source  5  in  a  self-absorbing  medium  of 
volume  Fand  density  p.  A  series  of  measurements  at  different  detector  positions  provide 
the  needed  information  on  the  source  distribution  and  its  total  activity.  Using  passive 
gamma  scanning  method,^  the  total  response  of  a  external  detector  at  position  /" ,  £)/ ,  to 
the  uncoUided  radiation  from  an  isotropic  source  is 

^'=h-s.-^^^-dv  (1.,) 

where  G/  is  a  local  constant  which  includes  all  geometrical  factors  and  the  detector 
efficiency  with  respect  to  detector  position  i;  Sy  is  the  volumetric  source  strength 
ir's/cm^  ■s);^is  the  linear  average  attenuation  coefficient;  tj  is  the  thickness  measured 
from  Sy  to  the  medium  boundary;  rj  is  the  distance  of  Sy  from  the  detector  at  position  /; 
and  dVis  the  volume  element.  Total  source  strength,  S,^,  (r's/s),  is  calculated  by 
integrating  Sy  over  the  container  volume. 

^tot  -  ]^y  -dV  (1-2) 

V 

The  inverse  radiation  transport  problem  which  is  considered  in  this  research  deals  with 
the  solution  of  a  series  of  integral  equations  similar  to  Equation  1-1  substituted  into 
Equation  1-2  for  the  calculation  of  total  source  strength  Sfot-  In  matrix  form,  the  problem 
can  be  represented  as  follows: 

D  =  KS  (1.3) 


Figure  1-1.  A  generic  configuration  of  a  distributed  source  S  in  a  self-absorbing  medium 
of  volume  V  measured  by  several  external  detectors 


where  D  is  the  vector  representing  the  external  detector  responses  (counting  rate),  K  is  the 
matrix  representing  the  point  response  kernel  and  all  the  geometrical  and  physical  factors 
involved  in  the  detection  process,  and  S  is  the  vector  representing  the  source  distribution 
in  the  drum.  In  most  cases,  the  distribution  of  the  radioactive  source  inside  a  nuclear 
waste  container  is  nonuniform.  For  a  near  field  measurement,  estimation  of  the  total 
activity  is  sensitive  to  (a)  spatial  resolution  (i.e.,  locations)  of  the  source,  (b) 
heterogeneity  of  the  absorbing  medium,  (c)  uncertainties  of  the  detecting  system,  and  (d) 
the  solution  scheme  used  to  determine  S.  Theoretically,  if  the  number  of  detector 
responses  is  equal  to  the  total  number  of  points  over  which  the  source  is  distributed, 
Equation  1-3  is  well  posed  and  an  exact  solution  can  be  obtained.  However,  due  to  the 
variability  of  the  source  distribution  and  the  heterogeneity  of  the  self -absorbing  medium. 
Equation  1-3  is  either  highly  singular  or  the  response  kernel  K  is  so  ill-conditioned  that  a 
unique  solution  can  not  be  found.  Uncertainties  involved  in  the  detection  process  and 
choice  of  solution  method  can  also  add  to  the  measurement  error.  Previous  attempts  to 
solve  this  problem  have  been  limited  to  cases  of  considering  a  single  point  source,  fixed 
source  distribution  patterns,  or  a  uniformly  distributed  source.^""  For  other  cases 
involving  a  wider  variety  of  source  distributions  and  density  heterogeneity,  a  satisfactory 
solution  to  this  problem  has  not  yet  been  found. 

This  research  develops  a  new  mathematical  approach  for  finding  an  optimum 
solution  to  a  singular  inverse  radiation  transport  problem  involving  the  nondestructive 
assay  of  radionuclide  inventory  in  a  large  nuclear  waste  container.  The  method  introduced 
is  an  optimizafion  scheme  based  on  performing  a  large  number  of  numerical  simulations 
which  account  for  the  counting  statistics,  the  nonuniformity  of  source  distribution  and  the 


heterogeneous  density  of  the  self-absorbing  medium  inside  the  nuclear  waste  container. 
The  simulation  model  uses  forward  projection  and  backward  reconstruction  algorithms. 
The  container  volume  is  first  modeled  to  have  several  computational  grids  for  a  vertical 
segment.  Then,  the  forward  projection  algorithm  uses  randomly  selected  source 
distribution  for  these  computational  grids  and  a  first  flight  kernel  method  to  calculate  the 
external  detector  responses.  The  backward  reconstruction  algorithm  uses  numerical 
optimization  methods  to  reconstruct  source  distribution  based  on  calculated  detector 
responses. '^"'^  Total  source  activity  is  determined  by  summing  the  reconstructed  activity 
of  each  computational  grid.  By  conducting  10,000  numerical  simulations,  the  error 
bound  and  the  associated  confidence  level  for  the  prediction  of  total  source  activity  is 
obtained.  A  relative  error  distribution  function  based  on  these  10,000  simulation  results 
is  also  determined.  The  accuracy  and  reliability  of  the  simulation  model  are  verified  by 
performing  a  series  of  experiments  in  a  208  I  waste  barrel.  A  number  of  Ir'^^  seeds  are 
used  to  simulate  unknown  distributions  of  the  source  in  the  barrel.  Density  heterogeneity 
is  simulated  by  using  different  materials  distributed  in  37  egg-crate  type  compartments 
simulating  a  vertical  segment  of  the  barrel.  Four  orthogonal  detector  positions  are  used  to 
measure  the  emerging  radiation  field  from  the  distributed  source. 

Chapter  2  describes  published  literature  related  to  the  nondestructive  assay  of 
radionuclide  inventory  in  a  nuclear  waste  drum.  Chapter  3  explains  the  modeling  of  the 
208  i  waste  barrel  and  the  idea  of  using  multiple  gamma  energy  peaks  for  the 
optimization  study  in  this  research.  Chapter  4  describes  the  conjugate  gradient  with 
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nonnegative  constraint  (CGNN)  reconstruction  method.'^  Chapter  5  describes  the 
maximum  likehhood  expectation  maximum  (MLEM)  reconstruction  method.'^ 

Chapter  6  details  experimental  setup,  procedures  of  data  collection,  and  results  of 
the  experiment.  Chapter  7  describes  concepts  and  procedures  for  performing  the 
numerical  simulations.  Results  of  the  activity  estimation  based  on  a  wide  spectrum  of 
simulated  experimental  conditions  are  presented.  Effects  of  the  density  heterogeneity  and 
locally  nonuniform  source  distribution  on  the  simulation  results  are  discussed.  A 
statistical  sampling  procedure  to  estimate  total  activity  for  the  full  length  of  the  waste 
barrel  is  demonstrated.  Chapter  8  presents  the  conclusions  of  this  research.  The  benefits 
and  the  potential  utilization  of  the  resuhs  of  this  research  are  also  presented. 


CHAPTER  2 
THE  NONDESTRUCTIVE  ASSAY  METHODS  FOR  CHARACTERIZING 
RADIONUCLIDE  INVENTORY  IN  A  NUCLEAR  WASTE  CONTAINER 

Segmented  Gamma  Scanning  Method 


The  segmented  gamma  scanning  (SGS)  method  was  developed  in  the  early 
1970s."  It  was  designed  to  assay  the  special  nuclear  waste  material  contained  in  low 
density  scrap  and  waste  by  measuring  gamma  emission  rates.  This  method  is  widely  used 
by  the  United  States  Department  of  Energy  complex  and  industry.  The  basic  assumption 
underlying  this  technique  is  that  materials  and  radionuclides  are  distributed 
homogeneously  and  uniformly  within  the  waste  barrel.  Barrels  are  usually  divided  into 
several  vertical  segments.  Gamma  measurement  is  done  one  at  a  time  for  each  vertical 
segment.  During  measurement,  the  waste  barrel  is  rotated  stepwisely  or  continuously  and 
the  uncollided  gamma  emission  rate  is  measured  by  a  collimated  germanium  detector 
located  outside  the  barrel.  In  a  continuous  scan,  the  rotating  sample  moves  past  the 
collimator  at  a  constant  speed.  In  a  stepwise  scan,  the  sample  is  positioned  vertically, 
counted,  repositioned,  counted  again,  and  so  on.  The  purpose  of  drum  rotation  is  to 
compensate  for  activity  distribution  asymmetries  such  that  the  uniform  and  homogeneous 
assumption  is  maintained.  Some  installations  use  an  external  gamma  source 
simultaneously  to  measure  the  average  density  and  correct  for  the  nonhomogeneous 
effect.  In  some  applications,  the  measurement  is  performed  at  a  very  large  distance 


(mostly  150  cm  to  300  cm  away  from  the  drum)  to  reduce  the  sensitivity  to  the  source 
distributions/"'  Such  applications  are  usually  limited  to  nuclear  waste  drums  with  much 
higher  activity  sources.  In  some  cases  involving  small  size  waste  barrels  or  low  activity 
levels,  detectors  are  positioned  at  close  distance  (within  30  cm)  and  homogeneous  source 
distribution  is  assumed.  For  the  estimation  of  total  activity  in  these  measurements,  the 
SGS  is  not  capable  of  making  an  accurate  assay  and  no  confidence  level  is  reported  for  a 
specific  high  accuracy  requirement. 

Tomographic  Gamma  Scanning  Method 

Since  1990,  many  efforts  have  been  devoted  to  the  tomographic  gamma  scanning 
(TGS)  method,  and  at  least,  two  groups  of  researchers  have  published  preliminary  results. 
One  was  a  Japanese  group  of  researchers  from  the  Hitachi  Ltd.  and  the  Toshiba 
Corporation.^"^  The  other  one  was  a  North  American  group  of  researchers  from  the 
Lawrence  Livermore  National  Laboratory  and  the  Los  Alamos  National  Laboratory.*"^ 
The  TGS  is  essentially  an  extension  of  the  SGS  method.  It  emphasizes  the  fact  that  the 
matrix  inside  the  drum  is  heterogeneous  in  nature,  and  the  density  distribution  of  the 
matrix  should  be  determined  by  a  transmission  tomographic  method  in  addition  to  the 
emission  spectrum.  A  representative  TGS  design  is  a  prototype  TGS  system  built  at  Los 
Alamos."^  It  uses  a  high  energy  gamma  transmission  source  moved  simuhaneously  with  a 
collimated  HPGe  detector.  For  the  horizontal  movement,  the  waste  drum  is  rotated 
stepwisely  or  continuously.  A  total  of  9  or  10  vertical  segment  scannings  are  made. 
Detailed  information  on  density  distribution  can  be  obtained  and  is  subsequently  fed  into 
the  traditional  SGS  method  to  correct  for  heterogeneous  matrix  effects. 


Compared  to  the  SGS  method,  TGS  takes  more  scamiing  and  analysis  time  in 
order  to  get  density  information.  For  nearly  homogeneous  waste,  sparsely  distributed 
solid  waste,  or  low  density  waste,  TGS  has  no  significant  advantage  over  SGS.  However, 
it  does  have  a  superior  advantage  if  matrix  density  is  heterogeneous  and  is  as  high  as 
2.0  g/cm^  High  density  waste  usually  consists  of  glassified  or  cemented  waste  forms  that 
are  nearly  homogeneous  in  density.  Therefore,  TGS  does  not  offer  much  improvement 
over  SGS,  unless  the  scanning  time  can  be  shortened  or  the  backward  reconstruction 
algorithm  can  be  improved. 

Passive  and  Active  Neutron  Detection  Method 

The  neutron  detection  method  is  specific  to  assays  of  waste  drums  containing 
actinides,  especially  for  the  measurement  of  plutonium  content.  Three  neutron  assay 
methods  are  reported:  (a)  the  passive  neutron  assay,  (b)  the  active  neutron  assay  using  an 
external  neutron  source,  and  (c)  the  active  neutron  assay  using  an  external  gamma 
source.'"""  The  passive  neutron  assay  method  utilizes  the  spontaneous  fission  behavior  of 
some  Pu  isotopes  (Pu^^^  Pu^^°,  and  Pu^^^)  and  detects  these  spontaneous  fission  neutrons 
by  He^  neutron  detectors.  The  active  neutron  assay  method  needs  an  external  source  to 
trigger  the  induced  fission  of  the  actinides.  External  sources  can  be  an  electron  linear 
accelerator  that  generates  high  energy  neutron  pulses  or  high  energy  gamma  rays,  or  it 
can  be  a  Cf "  spontaneous  fission  neutron  source.  When  using  external  neutron  pulses, 
these  fast  neutrons  are  thermalized  as  they  pass  through  a  moderator  matrix  surrounding 
the  waste  barrel  and  cause  fissile  actinides  to  produce  fast  fission  neutrons.  When  using 
external  gamma  rays,  high  energy  photons  (usually  as  high  as  7  MeV  or  above)  cause 
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production  of  fast  neutrons  from  all  actinides  by  photo-fission  and  photo-neutron 
reactions. 

Neutron  detection  methods  usually  need  a  coincidence  counting  system  and  time 
of  flight  detection  method  to  distinguish  between  time-correlated  fission  neutrons  emitted 
by  the  fission  process  and  random  neutrons  from  the  (a,n)  reaction  and  the  background. 
The  passive  neutron  assay  method  often  requires  a  long  measurement  time  (typically, 
several  hours  to  one  day  for  measuring  low  level  waste)  because  of  the  low  spontaneous 
fission  probability.  The  measurement  time  can  be  shortened  by  using  the  active  neutron 
assay,  but  the  accurate  determination  of  thermal  neutron  flux  within  the  moderator  matrix 
causes  another  measurement  difficulty.  Determination  of  the  inventory  of  actinides 
within  the  waste  drum  is  a  challenge  due  to  the  nature  of  nonhomogeneous  distribution  of 
actinides.  This  challenge  exists  whether  using  passive  or  active  neutron  assays  and  is  the 
same  challenge  faced  in  the  gamma  assay  method.  Therefore,  no  significant 
announcement  of  accurate  estimation  of  actinide  activity  inside  a  nuclear  waste  drum  has 
been  found  for  the  neutron  assay  method. 


CHAPTER  3 
THE  MULTIPLE  ENERGY  PEAKS  GAMMA  DETECTION  SYSTEM 

The  Model  of  a  Gamma  Scanning  System  Using  Four-Detector  Positions 

This  research  presents  a  calculation  model  that  uses  four  orthogonal  detector 
positions  to  measure  the  emerging  gamma  flux  from  the  distributed  source  in  a  vertical 
segment  of  a  208  ^  waste  barrel.  Figure  3-1  shows  the  configuration  of  a  vertical  segment 
of  the  waste  barrel.  Each  vertical  segment  is  modeled  to  have  37  computational  grids. 
Each  grid  is  assigned  a  number  as  indicated  in  Figure  3-1 .  The  volume  of  each  cubic  grid 
is  8x8x8  cm^  The  barrel  is  a  0.16  cm  thickness  steel  drum  and  has  56.84  cm  inner 
diameter.  The  gamma  signal  is  measured  at  a  distance  of  42.85  cm  from  the  center  of  the 
barrel.  The  calculation  model  assumes  that  the  HPGe  detector  (7.5  cm  in  diameter)  is 
well  collimated  and  detects  only  those  gamma  rays  coming  from  the  same  vertical 
segment. 

The  Concept  of  Using  Multiple  Gamma  Energv  Peaks 

Using  the  grid  system  as  described  above,  the  matrix  element  Kfj  in  Equation  1-3 
corresponds  to  the  counting  rate  of  a  specific  energy  peak  contributed  by  theyth  voxel  to 
the  detector  at  the  zth  position  and  is  given  by 

^/j-exp(-     X      {lL).p^.t^).cos{e)-s--^-p  (3-1) 

where  tq  is  the  path  length  of  the  ray(/j)  in  the  qth  voxel,  {m/p)  is  the  linear  mass 
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Figure  3-1 .  Configuration  of  a  four-detector  positions  measuring  system  in  a  vertical 
segment  of  a  208  ^nuclear  waste  barrel 
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attenuation  coefficient  of  the  gamma  energy  considered,  p^  is  ttie  density  in  the  ^th 
voxel,  Q,  ^  is  the  solid  angle  extended  from  a  source  point  in  theyth  voxel  to  the  detector 
at  the  rth  position,  s  is  the  detector  efficiency  at  a  calibration  position  with  solid  angle 
Qq  ,  and  /?is  the  branching  ratio  of  the  emitted  gamma  ray.  The  cos(6!)  in  Equation  3-1 
accounts  for  correction  of  the  inclination  if  the  gamma  ray  is  not  normal  to  the  detector 
surface.  In  real  implementation,  D  is  measured,  K  is  precalculated  based  on  a  known 
geometric  model  and  detector  efficiency,  and  S  can  be  solved  by  using  a  numerical 
optimization  algorithm  such  as  CGNN  or  MLEM.  Total  activity  can  be  obtained  by 
summing  over  all  the  elements  of  5*.  If  few  detector  positions  are  used,  this  is  an 
underdetermined  linear  algebra  problem.  If  single  peak  energy  is  used  to  analyze  this 
problem,  it  becomes  a  4x37  linear  system  for  which  a  reasonable  optimized  solution  is 
difficult  to  obtain.  Whenever  muhiple  energy  peaks  are  possible  for  analyzing  a  given 
radionuclide,  it  becomes  a  4mx37  linear  algebra  system  (m  being  the  number  of  identified 
energy  peaks).  If  the  number  of  detectors  is  increased  to  eight,  this  system  becomes 
8mx37.  For  the  four-detector  positions  and  two  energy  peaks  detection  system.  Equation 
1  -3  in  its  detailed  form  is 
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where  D,'"'^  is  the  mth  energy  peak  response  at  zth  detector  position,  J^/"''  is  the  point 
kernel  response  contributed  by  thejth  voxel  to  the  rth  detector  position,  and  5"^  is  the 
volumetric  source  in  theyth  voxel. 

This  research  uses  Ir"^  radionuclide  to  simulate  the  distributed  source.  Ir'^^  has 
several  prominent  and  well-separated  gamma  peaks  ranging  from  296  keV  to  884  keV 
(Figure  3-2),  and  is  appropriate  for  simulating  the  gamma  spectrum  of  plutonium  which 
originates  from  reprocessing  plants  and  is  of  great  interest  for  nuclear  waste  management. 
Three  energy  peaks  of  Ir''^  are  selected  for  this  research:  316.6  keV  (J3=  0.8285),  468.1 
keV  (J3=  0.4810),  and  588.6  keV  (p=  0.0457).  This  research  illustrates  that  from  a  two 
dimensional  analysis,  using  a  model  with  four-detector  positions  and  multiple  energy 
peaks  to  evaluate  the  activity  of  randomly  distributed  sources  inside  the  waste  drum, 
dramatic  improvement  of  the  accuracy  can  be  obtained.  In  the  numerical  simulation, 
source  strengths  are  randomly  generated  and  are  assigned  to  a  randomly  selected  grid.  If 
a  source  is  assigned  to  one  of  the  grids,  it  is  assumed  that  the  source  strength  is  uniformly 
distributed  within  that  grid. 

Computer  Modeling  and  Verification  Procedure 

The  following  steps  are  taken  for  the  computer  modeling  and  its  verification: 
1 .  A  computer  program  based  on  the  first  flight  kernel  method  is  written.  For 
every  measurement  or  numerical  simulation,  a  kernel  table  comprised  of  the  values  of  the 
elements  of  matrix  K  (shown  in  Equation  3-2)  is  generated.  Element  AT^J'  is  a  point 
kernel  response  which  is  an  averaged  value  of  the  counting  rate  (for  the  mth  energy  peak 
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Figure  3-2.  Ir     gamma  spectrum 
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at  the  rth  detector  position)  contributed  by  the  Monte  Carlo  sampling  of  5,000  point 
sources  within  grid;.  The  computer  program  can  generate  a  kernel  table  based  on  the 
consideration  of  heterogeneous  density  distribution. 

2.  Two  numerical  optimization  algorithms,  CGNN  and  MLEM,  are  built  into  the 
model  to  solve  the  linear  algebra  system.  The  convergence  criteria  are  selected  based  on 
the  study  of  the  optimum  result  from  10,000  random  source  distribution  cases. 

3.  The  computer  model  is  verified  by  various  physical  experiments  described  in 
Chapter  6.  These  experiments  involve  the  cases  of  nonuniform  distribution  of  Ir'^^ 
sources  under  heterogeneous  density  conditions. 

4.  Numerical  experiments,  described  in  detail  in  Chapter  7,  are  performed  to 


obtain  the  error  bound  and  the  associated  confidence  level  of  this  model  imder  the  i 

! 

influences  of  density  perturbation  and  signal  noise  perturbation.  First,  the  homogeneous 

condition  with  density  1 .0  g/cm^  (called  base  density)  is  analyzed  for  a  total  of  10,000  i 

r 
1 

I 

cases.  For  each  case,  total  source  activity  is  randomly  distributed  in  37  grids  and  each  I 

detector  response  for  each  energy  peak  is  calculated.  The  Monte  Carlo  procedure  is  used 

to  simulate  random  distribution  of  the  source  within  the  barrel.  A  random  number 

generator  is  used  to  randomly  choose  one  of  37  grids.  A  randomly  chosen  fraction  of 

total  source  activity  is  assigned  to  the  randomly  selected  grid.  Then,  a  second  grid  is 

chosen  and  a  randomly  selected  fraction  of  the  remaining  source  activity  is  assigned  to 

this  grid.  The  process  continues  until  all  remaining  source  activity  is  reduced  to  a  very 

small  value  (10"^  of  original  activity).  The  same  calculation  is  done  by  considering  the 

effect  of  density  heterogeneity.  Density  pertiorbation  is  simulated  by  selecting  a 

perturbation  range  which  is  some  fraction  (positive  or  negative)  deviated  from  the  base 
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density.  Again,  the  Monte  Carlo  procedure  is  used,  grid  by  grid,  to  randomly  choose  any 
value  within  the  range  and  assign  this  perturbed  density  for  the  grid.  The  detector 
responses  thus  calculated  are  called  the  noise-free  signal.  The  noise  perturbation  effect  is 
analyzed  by  considering  the  possible  uncertainties  of  detector  signal  in  a  real 
measurement,  which  may  be  caused  by  detector  calibration  uncertainties,  counting 
statistics,  and  the  locally  nonuniform  distribution  of  the  source  within  a  computational 
grid.  The  ideal  signal  is  perturbed,  first  by  a  selected  fraction  and  then  by  a  statistical 
sampling  process.  Several  noise  perturbation  levels  are  tested  and  analyzed. 


CHAPTER  4 
CONJUGATE  GRADIENT  WITH  NON-NEGATIVE  CONSTRAINT  METHOD 

The  conjugate  gradient  (CG)  algorithm  is  an  iterative  algebraic  solution  scheme 
for  solving  linear  systems  of  equations  and  represents  an  important  computational 
innovation  of  the  early  1950s.  This  method  came  into  wide  use  only  in  the  mid-1970s.''' 
Due  to  the  rapid  improvement  of  computer  technology,  iterative  algebra  algorithms, 
including  the  conjugate  gradient  method  and  expectation  maximum  (EM,  described  in  the 
next  chapter)  method,  have  become  more  and  more  important  in  the  fields  of  image 
reconstruction  and  remote  sensing  which  usually  require  a  lot  of  matrix  operations  and 
computer  memory.  An  application  similar  to  nuclear  waste  drum  characterization  is  the 
single  photon  emission  computed  tomography  (SPECT)  for  medical  imaging  in  clinical 
nuclear  medicine.  There  have  been  some  discussions  of  using  CG  in  SPECT,'^"'^  but 
wide  applications  have  not  been  seen.  SPECT  reconstructs  images  of  body  organs  which 
are  much  smaller  than  the  typical  nuclear  waste  drum.  The  total  activity  of  nuclear 
medicine  injected  into  the  human  body  is  always  known  and  is  usually  much  larger  than 
the  total  activity  in  waste  drum  characterization,  which  is  usually  unknown.  Therefore, 
applying  CG  in  the  nuclear  waste  drum  characterization  has  different  considerations  from 
that  in  the  SPECT. 

From  the  optimization  theorem,  the  least  squared  solution  S  of  Equation  1-3  can 
be  found  by  minimizing  the  following  squared  error  (see  the  derivation  in  Appendix  A): 
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D-KS 


=  {D-KS)\D-KS)  (4-1) 


where  |  |  means  the  L"  norm,  (  ^  means  the  transpose  of  the  matrix  and  S  means  the 

least  squared  solution  of  5.  In  the  conjugate  gradient  method,  solution  S  is  found  by 
using  an  iterative  algorithm: 

4+1  =  S„  +  a^d^  (4-2) 

where  5,,  is  the  estimate  of  S  at  the  nth  iteration,  a„  is  a  scalar  search  factor  and  J,,  is  a 
search  vector  which  is  required  to  be  orthogonal  to  the  modified  residual  vector  r„  by  the 
following  relation  : 

rLxd„  =  0  (4-3) 

with 

r„  =  K'{D-KS„)  (4-4) 

Equation  4-3  means  that  to  approach  the  minimum  point  of  e  ,  the  search  vector  at  a 
certain  iteration  step  must  be  orthogonal  to  the  residual  vector  at  the  next  step.  From  the 
relation 


^K'{KS„,,-D)  =  K\KS„~D)  +  a„K'Kd^ 
and  from  equation  (6),    we  get         r^^^  =  r^  -  a ^K' Kd ^  (4-5) 

To  minimize  e,  we  need  to  take  the  partial  derivative  of  Equation  4-1  with  respect 
to  a„  and  set  the  result  at  zero.  This  leads  to  (see  the  derivation  in  Appendix  B): 
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d{D-Ks")      c^(D-K{S„+a„dj')  ^'r 


The  search  vector  J„  of  CG  must  satisfy  the  general  Gram-Schmidt  conjugate  properties 
and  is  related  to  the  residual  vector  r  through  a  recursion  formula.'^    The  conjugate 
property  is 

d^K'Kd^  =0       ,      for  n^m  (4-7) 

and  the  recursion  relation  (see  the  derivation  in  Appendix  C)  is 

dn=r„-/3„^,d„_,    or<,,  =r„,i-^X  (4-8) 

where  /5„  is  the  orthogonality  factor. 

Taking  the  inner  product  of  Equation  4-8  with  K' Kd,,_^  and  applying  the  property 
of  Equation  4-7,  we  have 

The  algorithm  of  CG  without  constraint  is: 

Step  1 :  Set  an  initial  guess  of  5*, ,  and  calculate  d^^r^^  K'{D-KS^). 
Step  2:  Set  the  loop  index  starting  from  n=l,  calculate  «„  based  on 

Equation  4-6,  calculate  5„^i  based  on  Equation  4-2,  calculate  r„^i 
based  on  Equation  4-5,  calculate  /?„  based  on  Equation  4-9,  and 
then  calculate  <i„^,  based  on  Equation  4-8.  The  loop  ends  when  a 
specified  criterion  of  residual  difference  is  reached  or  a  specified 
iteration  number  is  reached. 
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The  algorithm  of  CGNN  is  more  complicated  than  the  above  algorithm.  Suppose 

there  is  a  hypersurface  //which  is  the  boundary  of  the  CGNN  solution  5,,^,  which 
minimizes  the  quadratic  function  e  of  Equation  4-1  and  has  at  least  m  components  which 
are  zero.  Then,  at  5*,,^, ,  any  slight  increase  of  component  5,™ ,  at  boundary  // will 
increase  the  e  value.  Therefore,  if  /  is  a  set  of  indices  i  <  m  such  that  5,',^i  ==  0 ,  we  have 
the  relations  for  the  optimum  point  5,,^,   : 

de  de 

^  U     for  i  in  /       and         -—^, ^    otherwise 


dS'  dS' 

Since  the  residual  is  defined  as  the  negative  value  of  the  gradient  (see  Equation  4-4  and 
Appendix  A),  the  criteria  for  terminating  the  search  and  finding  a  CGNN  solution  S*,,^;  is 
equivalent  to  have 

rj^i<0    for  fin/     and         ^^+i  =  0   otherwise  (4-10) 

The  search  method  is  conducted  by  a  restart  algorithm.  Within  any  iteration,  if  there 
exists  a  search  vector  dn  that  will  produce  negative  components  of  5'„^, ,  the  most 
negative  component  is  determined  and  the  scale  of  the  search  length  is  calculated  to  force 
this  component  to  be  zero.  Then,  the  /  set  is  redefined  and  the  criteria  in  Equation  4-10 

are  checked.  If  the  criteria  in  Equation  4-10  are  not  matched,   ^„^,  is  reset  to  be  S,  and 
the  whole  algorithm  is  restarted. 

The  most  important  issue  for  CGNN  algorithm  is  the  termination  point.  The  zero 
criterion  in  Equation  4-10  is  customarily  set  at  a  very  small  value  which  is  near  to  the 
precision  limit  of  the  computer  digit.  However,  because  the  source  activity  level  covers  a 
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very  wide  range,  a  fixed  value  for  zero  may  sometimes  stall  the  algorithm  and  may  also 
not  be  an  optimal  converged  point.  Zero  criterion  in  Equation  4-10  for  this  study  is 
dynamically  defined  to  be  a  small  fraction  of  the  average  value  of  the  initial  guess  of  r 
which  is  defined  in  Equation  4-4.  That  is 
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ZERO  =  ^ £  (4-11) 

37 

where  r,'"  is  the  mth  component  of  initial  guess  of  r  and  is  varied  case  by  case  with  respect 
to  the  measured  or  calculated  detector  responses.  Table  4-1  shows  the  convergence  test 
on  a  Pentium/120  computer  by  studying  10,000  random  source  distribution  cases  under 
homogeneous  density  condition.  A  total  activity  of  10*  Bq  is  input  for  every  case. 
Detector  responses  are  calculated  using  forward  projection.  Then,  the  source  distribution 
is  reconstructed  by  CGNN  algorithm  using  four  detector  positions  and  two  energy  peaks. 
The  reconstructed  total  activity  is  then  compared  with  the  input  total  activity  to  determine 
the  percent  error  for  each  case.    Five  categories  of  percent  error  are  defined  as  shovm  in 
Table  4-1.  The  cases  which  fall  into  a  certain  category  are  counted.  The  value  for  £-is 
varied  from  1 0"'  to  1 0"'°.  It  is  evident  that  the  optimal  value  of  £  among  them  is  1 0"*  which 
predicts  9,901  out  of  10,000  cases  with  10%  accuracy. 
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Table  4-1 .  Results  of  the  convergence  test  for  different  zero  criteria  of  CGNN  by  using 
10,000  random  source  distribution  cases  of  homogeneous  density 


e 

Error*  <50% 

<30% 

<25% 

<20% 

<10% 

10-' 

4,573 

2,957 

2,455 

2,057 

931 

10-^ 

8,356 

4,243 

3,418 

2,555 

1,207 

10-^ 

10,000 

9,929 

9,708 

9,213 

6,674 

10-' 

10,000 

10,000 

10,000 

9,993 

9,341 

10-^ 

10,000 

10,000 

10,000 

10,000 

9,783 

10-' 

10,000 

10,000 

10,000 

10,000 

9,901 

10-' 

10,000 

10,000 

10,000 

10,000 

9,823 

10-^ 

10,000 

10,000 

10,000 

10,000 

9,824 

10-' 

10,000 

10,000 

10,000 

10,000 

9,823 

10-'° 

10,000 

10,000 

10,000 

10,000 

9,820 

*  Error  means  the  percent  estimation  error  of  the  total  activity  and 
the  values  in  the  2nd  to  5th  column  are  the  number  of  cases  which 
fall  into  that  error  category. 


CHAPTER  5 
MAXIMUM  LIKELIHOOD  EXPECTATION  MAXIMUM  METHOD 

Maximum  likelihood  (ML)  method  is  an  approach  to  maximize  the  likelihood 
function,  which  is  the  probability  that  the  source  strength  S  produces  the  measured 
detector  responses.  As  is  CG,  ML  is  an  iterative  algebraic  solution  scheme  for  solving 
linear  systems  of  equations.  It  was  introduced  to  emission  tomography  by  Rockmore  and 
Macovski  in  1976."  Lange  and  Carson  incorporated  the  expectation  maximum  (EM) 
method  in  ML  to  compute  maximum  likelihood  estimates."  MLEM  then  became  a 
popular  application  on  SPECT.  It  assumes  that  both  £>  and  -S  in  Equation  1-3  are  Poisson 
distributions.  Thus,  the  expectation  value  of  the  ith  detector  response  for  a  waste  barrel 
with  37  computational  grids  is 

D.=tK.A=t^,  (5-1) 

where  5^-  is  the  source  strength  at  voxel;  and  Ky  is  the  probability  that  a  photon  leaving 
voxel  J  is  counted  by  the  /th  detector.  X,^  is  the  expectation  value  of  Xy,  which  is  the 
statistical  number  of  photons  that  are  emitted  from  voxel/  and  contributes  to  the  counts 
for  detector  response  /.  In  the  EM  algorithm,  a  conditional  expectation  function  is 
defined  as 

E(\nf{X,S)\D,S^"^)  (5-2) 

where  In  denotes  the  natural  logarithm  and  5"^"'  is  the  estimation  of  S  at  the  nth  iteration. 
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This  equation  defines  an  expectation  function  of  a  logarithm  likelihood  function /(X,^) 
under  the  condition  of  observed  detector  response  D  and  source  distribution  S  at  the  nth 
iteration  step.  f{X,S)  has  the  Poisson  distribution  form  of 

m      37     X^"   -e'^'' 
/=1    ;=1  ^  ij  • 

where  m  is  the  total  number  of  detector  response.  Its  logarithm  is 
\nf{X,S)  =  22:(-^./  +  ^.  1^^.  -  ^^^. ') 

=  Z  S  (-^.'^. + ^ij  ^<^u^j )  - 1^  ^^ ')  (5-4) 


'        7 


The  conditional  expectation  of  X^  in  Equation  5-4  with  respective  to  £)/  and  5^"'  is 

K.S^"^  ■  D. 
N,=EiX,\D^,S^"^)--^^^  (5-5) 


Because  the  third  term  in  Equation  5-4  has  no  dependence  on  the  iteration  of  5,  the 
conditional  expectation  of  Equation  5-4  becomes 

E(\nfiX,S)\D,S'"')  =  Ut-^.^.  +  ^.  HK,Sj)]  +  C  (5-6) 

where  C  is  the  term  not  dependent  on  the  iteration  ofS.  To  find  the  maximum  of  this 
expectation,  we  need 

^    Eilnf(X,S)\D,S^"^)  ^  ^         1 


=  0  =  -S^,+E^.^  (5-7) 


^    S^  ^    '^     r     "Sj 


The  partial  derivative  of  the  above  equation  is 
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d    '     E{\nf{X,S)D,S("')  1 

=  -I^.7;7  (5-8) 


d    s;  ^    '  S] 


Equation  5-8  yields  a  negative  value  and  this  assures  that  E{\nf{X,S)  D,S  "  )  is  a 
concave  function  of  S  and  has  a  maximum  at 

'    ~  Tk,  ~  JlK,  f  2;x„sr'   •  y,/=i....,37 

This  iteration  algorithm  will  always  assure  nonnegative  results  if  the  initial  guess  of  the 
source  is  positive. 

To  determine  the  optimal  termination  point  for  the  iterative  algorithm  in  Equation 
5-9,  the  same  procedvire  as  described  in  the  last  paragraph  of  Chapter  4  is  adopted.  The 
termination  point  is  set  to  be  the  optimal  number  of  iterations  which  will  give  the  most 
satisfactory  results  from  the  process  of  studying  10,000  random  source  distribution  cases. 
Table  5-1  shows  the  results  of  the  convergence  test  by  varying  the  iteration  number  from 
100  to  20,000.  The  optimal  prediction  among  them  is  10,000  iterations.  MLEM  predicts 
9,350  out  of  10,000  cases  with  10%  accuracy.  The  computation  time  of  MLEM  using 
10,000  iterations  is  about  10  times  longer  than  that  of  CGNN  using  zero  criterion  10"*. 
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Table  5-1.  Results  of  the  convergence  test  for  different  iteration  number  of  MLEM  by 
using  10,000  random  source  distribution  cases  of  homogeneous  density 


It.#     Error*  <50% 

<30% 

<25% 

<20% 

<10% 

100 

9,906 

6,918 

5,832 

4,572 

2,821 

500 

10,000 

9,963 

9,872 

9,582 

6,547 

1,000 

10,000 

10,000 

9,991 

9,831 

7,800 

2,000 

10,000 

10,000 

10,000 

9,906 

8,573 

5,000 

10,000 

10,000 

10,000 

9,942 

9,122 

10,000 

10,000 

10,000 

10,000 

9,943 

9,350 

20,000 

10,000 

10,000 

10,000 

9,943 

9,163 

*  Error  means  the  percent  estimation  error  of  the  total  activity  and 
the  values  in  the  2nd  to  5th  column  are  the  number  of  cases  which 
fall  into  that  error  category. 
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CHAPTER  6 
EXPERIMENTS  AND  MODEL  VERIFICATIONS 

A  series  of  experiments  were  performed  to  verify  the  computer  model  and  to 
evaluate  the  properties  of  the  concept  of  using  multiple  energy  peaks  in  numerical 
reconstruction  algorithms.  These  experiments  were  designed  to  simulate  the  random 
source  distributions  and  the  heterogeneous  density  conditions  which  are  usually  met  in 
real  waste  barrel  measurements.  To  simulate  the  distributed  source,  several  pieces  of  Ir 
sources  were  prepared.  Each  piece  of  Ir'^^  source  was  3  mm  long  and  0.5  mm  in 
diameter  and  had  activity  of  0.0726  ±  5%  mCi  (2.6869x10*  ±  5%  Bq)  on  October  16, 
1996.  Three  Regions  of  Interest  (ROI)  in  the  Ir'^^  gamma  spectrum  were  defined  for  the 
peak  area  calculation.  These  ROIs  were  centered  at  3 16.6  keV  iJ3=  0.8285),  468.1  keV 
{J3=  0.4810),  and  588.6  keV  (y9=  0.0457)  respectively.  For  reconstruction  using  one 
energy  peak,  the  3 16.6  keV  peak  was  used.  For  reconstruction  using  two  energy  peaks, 
the  3 16.6  keV  and  588.6  keV  peaks  were  used;  and  for  reconstruction  using  three  energy 
peaks,  all  three  peaks  were  used. 

Experimental  Setup 

Figure  6-1  is  a  photograph  of  the  experimental  setup  including  a  208  i  waste 
barrel,  lead  collimators,  a  HPGe  gamma  detector,  electronic  modules,  and  a  data  access 
computer.  The  HPGe  gamma  detector  (EG&G  ORTEC  model  GEM-10195-P)  has  a 
transistor-reset  type  preamplifier  and  has  10%  relative  efficiency  and  1.95  keV  resolution 
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measured  at  1 .33  MeV.  The  electronic  modules  include  an  amplifier  (EG&G  ORTEC 
model  572),  a  bias  voltage  supply  (EG&G  ORTEC  model  459),  and  a  multichannel 
analyzer  (MCA)  interface  card  (EG&G  ORTEC  920  MCB). 

The  waste  barrel  had  a  steel  wall  0.16  cm  thick.  Its  inner  diameter  was  measured 
56.84  cm  and  the  height  was  88  cm.  The  waste  barrel  was  seated  on  top  of  a  rotation 
device  which  was  capable  of  360  degree  manual  rotation.  A  cylindrical  wooden  stand 
25.4  cm  high  and  with  about  the  same  inner  diameter  as  the  barrel  was  built  and  placed 
inside  and  at  the  bottom  of  the  barrel.  To  simulate  the  37  computational  grids  for  a 
vertical  segment  of  the  barrel,  a  matrix  container  (see  Figure  6-2)  with  37  egg-crate  type 
compartments  made  of  acrylic-glass  (1.5  mm  in  thickness)  was  built  and  cut  to  fit  the 
round  shape  of  the  barrel  at  the  periphery.  This  acrylic-glass  matrix  container  was  seated 
on  top  of  the  wooden  stand.  Each  of  the  37  compartments  was  8x8  cm^  in  area  and  25.4 
cm  high.  The  HPGe  detector  was  located  42.85  cm  away  from  the  center  of  the  barrel 
(14.31  cm  away  firom  the  periphery  of  the  barrel)  with  its  center  line  at  half  of  the  height 
of  the  matrix  container,  and  coUimated  by  three  lead  bricks.  The  firont  brick  had  an  open 
hole  about  the  same  diameter  as  the  face  of  the  detector  and  allowed  the  detector  to 
penetrate  through  it.  Two  other  lead  bricks  were  positioned  behind  the  front  one  and 
seated  on  both  sides  of  the  detector  to  ensure  coverage  of  the  whole  collimating  system. 

The  data  access  computer  was  a  PC-AT  with  a  MCA  interface  card  installed. 
The  gamma  ray  signal  was  processed  by  electronic  modules  and  signal  counts  were 
stored  on  the  buffer  of  the  MCA  card.  EG&G  ORTEC  software  MAESTRO-II  controlled 
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Figure  6-1.  Experimental  setup  for  the  208  ^  waste  barrel  characterization 
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Figure  6-2.  The  acrylic-glass  matrix  container  inside  the  waste  barrel 
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the  timing  of  the  signal  processing  and  interpreted  the  information  stored  on  the  MCA 
buffer  as  pulse  height  data  (often  called  an  energy  spectrum). 

Detector  Efficiency  Calibration 

Five  check  sources  were  used  for  HPGe  detector  efficiency  calibration.  The 
source  activities  on  September  15,  1996,  were  1.037  |aCi,  1.025  |aCi,  0.9173  [iCi,  0.9084 
^Ci,  and  0.8978  ^iCi  for  Ba'",  Cs'",  Mn'',  Na'',  and  Co"  respectively.  Each  check 
source  was  placed  along  the  center  line  of  the  detector  and  at  38.1  cm  (15  inches)  away 
from  the  detector  front  face.  Electronic  settings  for  the  calibration  were  bias  voltage 
2300V,  amplifier  gain  24,  and  shaping  time  2  )j.sec.  Table  6-1  shows  the  characteristics 

i 

of  these  five  check  sources  and  the  measured  gamma  energy  efficiencies.  The  energy  \ 

i 
I 

range  extended  from  8 1 .2  kev  to  1 33 1 .0  kev.  Figure  6-3  shows  the  measured  efficiencies  \ 

and  the  resulting  calibration  curve  based  on  least  squares  fitting.  Two  steps  of  least  ; 

squares  fitting  were  used  to  establish  the  calibration  curve.  A  second  order  polynomial  { 

\ 
r 

fitting  was  used  from  8 1 .2  to  136.8  keV  and  a  third  order  polynomial  fitting  was  used  [ 

[ 
from  136.8  to  1331  keV.  The  equation  for  the  third  order  polynomial  fitted  curve  is  i 

e(E)=  4.44942  x  10^'  -  9.88770  x  10"'  £  +  9.60628  x  lO"'"  £^  -  3.22253  x  10"'^  £'  1 

where  s  is  the  point  source  efficiency  at  a  distance  of  1 5"  from  the  detector,  and  E  is  the  \ 

gamma  energy  in  keV.  Efficiencies  based  on  this  calibration  curve  for  3 16.68  keV,  ' 

468.09  keV,  and  588.09  keV  energy  peaks  of  Ir'''  source,  were  found  to  be  2.17921x10"^,  f 

1.59539x10"^,  and  1.30146x10"^  respectively.  \ 
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Table  6-1 .  Characteristics  of  check  sources  and  related  efficiency  calibrations 


Gamma  Energy  for  the  Efficiency  Calibrated  at 

Check  Sources  (key)        Nuclide*  (Branching  Ratio)      15"  from  the  Detector 


81.2 
122.3 
136.8 
276.3 
302.8 
356.0 
383.9 
662.0 
835.1 
1173.0 
1331.0 


Ba'" 

:0.3297) 

2.1510x10-' 

Co" 

:0.8560) 

3.4730x10-" 

Co" 

:0.1068) 

3.2070x10-" 

Ba'^^  < 

:0.0715) 

2.7140x10-" 

Ba'^^  ( 

:0.1830) 

2.2530x10-" 

Ba'^^  ( 

;0.6194) 

1.8405x10-" 

Ba'^^  ( 

;0.0891) 

1.6799x10-" 

Cs'"  ( 

0.8510) 

1.2885x10-" 

Mn'*  ( 

;0.9998) 

1.0325x10-" 

Co*"    ( 

;i.oooo) 

7.9700x10-' 

Co*"    ( 

;i.oooo) 

7.4400x10-' 

*  Half-life:  Ba"' -10.57  years,  Cs'"-30. 17  years,  Mn'"-312  days, 
Na^^-2.60  years,  Co"-271.8  days. 
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Figure  6-3.  The  efficiency  calibration  curve  for  HPGe  detector 
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Region  of  Interest  of  Ir''^^  Gamma  Spectrum 

MAESTRO-II  calculates  total  counts  of  a  measured  energy  peak  by  integrating 
over  the  peak  area  defined  by  ROI.  ^°  Precise  definition  of  the  ROI  is  very  important  to 
accurate  peak  area  calculation.  Figure  6-4  shows  the  configuration  for  defining  a  ROI  of 
a  316  kev  energy  peak.  The  ROI  is  a  group  of  MCA  channels  bounded  by  a  low  channel 
number  /  and  a  high  channel  number  h.  The  background  on  the  low  channel  side  of  the 
peak  is  the  average  of  the  first  three  channels  of  the  ROI.  The  background  on  the  high 
charmel  side  of  the  peak  is  the  average  of  the  last  three  channels  of  the  ROI.  These  two 
pairs  of  channels  define  the  end  points  of  the  straight-line  background.  The  background 
is  given  by  the  following  equation: 

V,W  ,=A-2      ^  6 

where  B  is  the  total  background  area;  /  is  the  ROI  low  limit;  h  is  the  ROI  high  limit;  and 
Q  is  the  number  of  counts  of  channel  /.  A  gross  area  Ag  is  defined  as 

4=EC,  (6-2). 

/=/ 

The  adjusted  gross  area  is  the  sum  of  all  the  channels  within  the  ROI  but  not  used  in  the 

background,  and  is  represented  by  the  following  relation: 

A, -EQ  (6-3) 

The  net  area  is  the  adjusted  gross  area  minus  the  adjusted  background  as  follows: 
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Figure  6-4.  The  ROI  for  a  3 16  keV  energy  peak  of  the  Ir"^  gamma  spectrum 
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This  research  uses  the  net  area  counts  A^  as  the  total  uncolhded  photons  measured  by  the 
detector  for  a  specific  gamma  energy  peak. 

The  ROIs  of  Ir'^-  were  defined  at  the  time  that  the  energies  were  calibrated.  As 
shown  in  Figure  6-4,  the  energy  calibration  set  the  center  of  316  keV  peak  at  channel  775, 
with  a  lower  ROI  limit  set  at  channel  766  (3 1 3.52  keV)  and  a  higher  ROI  limit  set  at 
channel  782  (319.13  keV).  The  Full- Width-Half-Maximum  (FWHM)  of  this  peak  is  1.77 
keV.  The  ROI  encompasses  a  region  three  times  that  of  the  FWHM.  Figure  6-5  shows 
the  configurations  of  the  ROI  for  the  468  keV  and  588  keV  energy  peaks  based  on  the 
calibrated  energy  spectrum.  The  468  keV  energy  peak  was  centered  at  channel  1 150  and 
had  a  FWHM  of  1.85  keV,  a  low  ROI  limit  of  1141,  and  a  high  ROI  limit  of  1 157.  The 
588  keV  energy  peak  was  centered  at  channel  1448  and  had  a  FWHM  of  1 .87  kev,  a  low 
ROI  limit  of  1440,  and  a  high  ROI  limit  of  1456. 

Data  Collection  Process 

To  simulate  a  four  fixed  detector  system,  detector  readings  were  taken  at  each  90 
degree  position  by  rotating  the  barrel  counterclockwise.  For  the  homogeneous  condition, 
the  barrel  was  flooded  with  water  and  18  nonuniform  source  distribution  cases  were 
studied.  For  heterogeneous  conditions,  in  addition  to  water,  six  other  materials  were 
used:  sand  (p  =  1.45  g/cm'),  small  stones  (p  =  1.33  g/cm^),  sah  (p  =  1.14  g/cm^),  flour  (p 
=  0.83  g/cm^),  bird  seed  (p  =  0.78  g/cm^),  and  pine  ashes  (p  =  0.27  g/cm^).  Figure  6-6 
shows  that  the  acrylic-glass  insertion  chaimels  were  specifically  designed  to  hold  the 
above  mentioned  materials  and  to  exactly  fit  into  the  37  compartments.  A  thin  string  was 
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Figure  6-5.  The  ROI  for  the  468  keV  and  588  keV  energy  peaks  of  the  Ir"^  spectrum 


attached  to  the  top  of  every  insertion  channel  to  facilitate  its  insertion  into  and  removal 
from  the  compartment.  Figure  6-6  also  shows  a  configuration  of  an  experimental 
heterogeneous  condition  inside  the  waste  barrel.  A  number  of  insertion  channels  are  used 
in  this  configuration.  Plastic  bags  were  used  for  water  channels  to  avoid  the  possibility  of 
leakage.  Different  density  distributions  were  modeled  by  shuffling  the  insertion 
chaimels.  Two  heterogeneous  conditions  (Type  A  and  Type  B)  were  simulated  by  using 
combinations  of  the  seven  materials.  Density  Type  A  simulated  low  density  perturbation 
condition  and  its  density  distribution  varied  from  0.83  to  1.14  g/cm".  Density  Type  B 
simulated  high  density  perturbation  condition  and  its  density  distribution  varied  from  I 

0.27  to  1 .45  g/cml  Both  Type  A  and  Type  B  conditions  had  9  nonuniform  source  j, 

! 

distribution  cases. 

For  each  source  distribution  case,  Ir"^  sources  were  placed  at  the  center  of  the 
height  of  some  insertion  channels.  In  order  to  avoid  any  bias  resulting  from  orientation 
of  Ir'^^  sources,  and  to  obtain  reasonable  statistical  results,  four  300  second  measurements 
were  made  at  each  detector  position.  Insertion  channels  containing  Ir'^^  sources  were 
removed,  rotated  90  degree  clockwise,  and  reinserted  into  their  compartments  for  each 
successive  measurements.  Detector  responses  were  then  averaged  for  these  four 
measurements. 

Experiments  of  Distributed  Ir'^^  Sources  under  the  Homogeneous  Density  Condition 


Figure  6-7  shows  the  positions  of  Ir"'^  sources  for  18  experiments  under  the 
homogeneous  density  condition.  These  positions  were  chosen  to  represent  different 


40 


Figure  6-6.  Six  materials  of  different  density  in  the  insertion  channels  and  an 
experimental  heterogeneous  configuration  inside  the  waste  barrel 
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source  distribution  behaviors.  Cases  HI,  H2,  H3,  and  HI 3  represented  single  isolated 
source  distribution.  Cases  H7,  H8,  H9,  HI  1,  H14,  H15,  H16,  and  H17  represented 
clustered  source  distributions.  All  other  cases  belonged  to  multiple  but  scattered  source 
distributions.  In  the  figure,  the  activity  marked  on  top  of  each  configuration  was  the  total 
activity  for  that  configuration.  Case  HI 4  had  three  Ir''^  pieces  located  at  grid  20  and 
another  piece  at  grid  21.  Case  HI  5  had  three  Ir'^^  pieces  located  at  grid  25  and  three  other 
pieces  scattered  at  grid  6,  14,  and  33  respectively. 

Four  measurements  were  made  at  each  detector  position.  For  each  measurement, 
total  counts  of  316  keV,  468  keV,  and  588  keV  energy  peaks  were  recorded  respectively, 
and  then  these  values  were  divided  by  counting  time  to  get  detector  responses  (counting 
rates).  Final  detector  responses  were  averaged  from  these  four  measurement  values. 
Four  orthogonal  detector  positions  were  measured  with  a  total  of  16  measurements  made 
for  each  experiment.  These  measured  detector  responses  formed  the  elements  of 
response  matrix  D  in  Equation  1-3,  and  when  combined  with  the  precalculated  kernel 
matrix  K,  a  numerical  reconstruction  algorithm  was  used  to  solve  for  source  distribution. 
When  two  energies  were  used  for  the  reconstruction,  the  linear  system  was  8x37.  When 
three  energies  were  used,  the  linear  system  was  12x37.  Results  of  total  activity 
estimation  for  these  18  experiments  are  shovm  in  Table  6-2.  Using  the  CGNN  method, 
estimation  errors  ranged  from  -33.30%  to  27.52%  for  2-energy  reconstruction  and  from 
-34.38%)  to  28.67%  for  3-energy  reconstruction.  Using  the  MLEM  method,  estimation 
errors  ranged  from  -29.22%  to  23.79%  for  2-energy  reconstruction  and  from  -28.30%  to 
24.62%  for  3-energy  reconstruction.  Table  6-3  shows  the  summary  of  these  results.  All 
of  the  18  experiments  were  within  50%  error  for  the  four  reconstruction  schemes  used. 
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Figure  6-7.  Configurations  of  Ir'^^  source  distribution  for  18  experiments  under  the 
homogeneous  density  condition  (®  :  1  Ir''-  piece  ;  •  :  3  Ir'^^  pieces) 
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Table  6-2.  Activity  estimation  errors  for  18  homogeneous  density  experiments  by  using 
multiple  energy  peaks  in  the  CGNN  and  MLEM  reconstruction  algorithms 


CGNN 

MLEM 

Reconstruction 

Reconstruction 

Number  of 

Error  (%) 

Error  (%) 

Ir"'^  Pieces 

Case 

2 -Energy 

3 -Energy 

2 -Energy 

3 -Energy 

in  the  Barrel 

HI 

-3.90 

-2.09 

-3.98 

-1.85 

1 

H2 

-9.54 

-8.38 

-9.65 

-8.55 

1 

H3 

27.52 

28.67 

-1.37 

-1.40 

1 

H4 

4.91 

2.84 

5.11 

6.56 

2 

H5 

4.21 

3.52 

5.21 

4.13 

2 

H6 

-14.29 

-11.76 

-14.92 

-12.29 

2 

H7 

-1.14 

-0.16 

-2.40 

-1.03 

2 

H8 

3.24 

3.44 

1.36 

2.49 

3 

H9 

15.07 

13.92 

-1.26 

-0.61 

3 

HIO 

21.32 

23.65 

23.79 

24.62 

3 

Hll 

3.58 

4.55 

1.72 

3.53 

4 

H12 

-6.58 

-6.35 

-8.41 

-7.86 

2 

H13 

-21.10 

-18.51 

-20.59 

-15.58 

1 

H14 

-25.50 

-34.38 

-28.15 

-25.52 

4 

H15 

-33.49 

-31.55 

-29.22 

-28.30 

6 

H16 

-30.63 

-29.70 

-26.55 

-25.44 

4 

H17 

-17.56 

-14.17 

-19.12 

-18.42 

4 

H18 

-29.51 

-29.39 

-26.43 

-24.77 

4 

44 


Table  6-3.  Summary  of  activity  estimation  errors  for  18  homogeneous  density 
experiments  by  using  multiple  energy  peaks  in  the  CGNN  and  MLEM  algorithms 


Error  <50%         <30%         <25%         <20%         <10% 


2-Energv  Reconstruction 


CGNN          18'             17               13  11  8 

(100%)**  (94.4%)  (76.5%)  (61.1%)  (44.4%) 

MLEM          18              17               13  12  8 

(100%)  (94.4%)  (76.5%)  (66.7%)  (44.4%) 

3 -Energy  Reconstruction 

CGNN          18              18               13  12  8 

(100%)  (100%)  (76.5%)  (66.7%)  (44.4%) 

MLEM          18              18               15  13  10 

(100%)  (100%)  (83.3%)  (76.5%)  (55.6%) 


Number  of  cases  within  the  specified  error  bound. 
Percentage  of  the  1 8  cases  within  the  specified  error  bound. 


45 

For  30%  error,  both  the  CGNN  and  MLEM  predicted  17  out  of  18  cases  (94.4%)  by  using 
2-energy  reconstruction,  and  18  out  of  18  cases  (100%)  by  using  3 -energy  reconstruction. 
For  20%  error,  the  best  prediction  of  the  four  reconstruction  schemes  was  the  MLEM 
using  3-energy  reconstruction  which  predicted  13  out  of  18  cases  (76.5%). 

Experiments  of  Distributed  Ir'^^  Sources  under  Heterogeneous  Density  Conditions 

Two  types  of  density  heterogeneity  were  studied.  Type  A  heterogeneity 
simulated  a  low  density  perturbation  condition  by  filling  insertion  channels  with  salt, 
water,  or  flour  and  then  inserting  them  into  the  egg-crate  compartments  in  the  waste 
barrel.  Type  B  heterogeneity  simulated  a  high  density  perturbation  condition  by  using 
sand,  small  stones,  salt,  water,  flour,  bird  seed  and  pine  ashes.  Before  placing  these 
materials  into  the  barrel,  their  densities  were  measured  by  transmission  density 
measurement.  An  Ir"-  piece  was  used  as  the  transmission  source  and  was  placed  15" 
away  from  the  detector.  Counting  rates  of  588  keV  peak  were  recorded  for  measurements 
both  with  and  without  an  insertion  charmel  in  front  of  the  detector.  Density  then  was 
calculated  based  on  the  exponential  attenuation  relation.  For  each  material,  a 
transmission  density  measurement  was  made  for  each  of  the  four  faces  of  the  rectangular 
charmel.  The  density  for  each  material  listed  in  Table  6-4  was  the  average  value  from 
these  four  measurements.  Table  6-4  also  lists  the  positions  of  different  material  channels 
inserted  into  the  barrel  for  both  Type  A  and  Type  B  heterogeneity  conditions.  Figure  6-8 
shows  the  configuration  of  Ir'^^  locations  for  the  nine  experiments  under  the  Type  A 
heterogeneity  condition.  Source  configurations  for  the  other  nine  experiments  under  the 
Type  B  heterogeneity  condition  are  shown  in  Figure  6-9.  As  in  the  homogeneous  density 
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Table  6-4.  Location  of  different  material  channels  inside  the  barrel  for  the  study  of  Type 
A  and  Type  B  density  heterogeneity 


Type  A  Heterogeneity  (Density  :  0.83  -1.14  g/cm^) 


Grid  Number 


Salt    (p=1.14g/cm^)  17,19,21,25,27,28,31,33,34,36 

Water  (p  =1.00  g/cm')  1,    2,    3,    9,12,14,15,16,18,20, 

22,  23,  26,  29,  32,  35,  37 
Flour  (p  =0.83  g/cm')  4,    5,    6,    7,    8,10,11,13,24,30 


Type  B  Heterogeneity  (Density  :  0.27  ~  1.45  g/cm^) 


Grid  Number 

Sand  (p  =1.45  g/cm')  2,  4,  5,  6,  1 1 

Stones  (p  =1.33  g/cm^)  10,  17,  18 

Salt  (p  =1 .14  g/cm')  8,  25,  26,  3 1,  32,  36 

Water  (p  =1.00  g/cm')  1,3,    9,15,16,22,23,29,35,37 

Flour  (p  =0.83  g/cm')  7,  12,  13,  14,  20,  33,  34 

Bird  Seed  (p  =0.78  g/cm^)  19,24,30 

Pine  Ashes  (p  =0.27  g/cm^)  21,  27,  28 
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Figure  6-8.  Density  distribution  and  configurations  of  Ir''^  location  for  9  experiments 
under  the  Type  A  heterogeneity  condition 
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Figure  6-9.  Density  distribution  and  configurations  of  Ir'^^  location  for  9  experiments 
under  the  Type  B  heterogeneity  condition 
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cases,  the  Ir'^^  source  locations  (shown  in  these  two  figure)  simulate  single,  clustered, 
and/or  multiple  but  scattered  source  distributions. 

The  resuhs  of  total  activity  estimation  for  these  experiments  appear  in  Table  6-5. 
For  all  of  the  reconstruction  schemes  used,  the  activity  of  most  of  the  cases  is 
underestimated.  The  underestimation  may  be  attributed  to  the  biases  of  the  average 
density  measurements  of  the  insertion  channels,  or  to  the  biases  of  the  choice  of  Ir'^^ 
distribution.  Since  the  number  of  experiments  is  limited,  the  exact  bias  tendency  of  the 
estimation  can  be  assured  only  after  performing  more  experiments  through  the  numerical 
simulations. 

Table  6-6  summarizes  experimental  results.  All  of  the  experiments  were  within 
50%  error  for  the  four  reconstruction  schemes  used.  In  Type  A  heterogeneity 
experiments  with  30%  error,  CGNN  predicted  7  out  of  9  cases  (77.8%)  by  using 
2-energy  reconstruction,  and  MLEM  predicted  8  out  of  9  cases  (88.9%).  By  using 
3-energy  reconstruction,  both  CGNN  and  MLEM  predicted  8  out  of  9  cases  (88.9%).  In 
Type  B  heterogeneity  experiments  with  30%  error,  CGNN  predicted  8  out  of  9  cases 
(88.9%)  by  using  2-energy  or  3-energy  reconstruction,  and  MLEM  predicted  9  out  of  9 
cases  (100%)  by  using  2-energy  or  3-energy  reconstruction.  All  the  9  cases  predicted  by 
MLEM  using  3-energy  reconstruction  were  within  20%  error. 

Discussion  of  the  Results  of  the  Experiments 

Two  important  findings  were  observed  from  the  above  36  experiments:  (a)  the 
computer  model  was  able  to  accurately  estimate  total  activity,  and  (b)  the  locally 
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Table  6-5.  Activity  estimation  errors  for  the  heterogeneous  density  experiments  by  using 
multiple  energy  peaks  in  the  CGNN  and  MLEM  algorithms 


CGNN 

MLEM 

Reconstruction 

Reconstruction 

Number  of 

Error  (%) 

Erroi 

•  (%) 

Ir'^' Pieces 

Case 

2 -Energy 

3 -Energy               2 

-Energy 

3 -Energy 

in  the  Barrel 

Cases  ot  lype  A  Heterogeneity:  (A1-A9) 

Al 

-38.77 

-34.18 

-31.66 

-30.15 

2 

A2 

-18.36 

-15.34 

-19.56 

-16.44 

1 

A3 

-36.86 

-16.61 

-23.74 

-12.32 

2 

A4 

1.36 

4.49 

2.16 

5.78 

1 

A5 

3.74 

3.13 

0.95 

1.23 

1 

A6 

-23.75 

-14.72 

-25.09 

-17.23 

2 

A7 

-24.42 

-22.12 

-19.99 

-18.25 

3 

A8 

-9.40 

-12.67 

1.96 

1.55 

3 

A9 

-19.01 

-20.66 

-17.34 

-15.38 

4 

Cases  of 

Type  B  Heterogeneity:  (B1-B9) 

Bl 

-24.74 

-22.82 

-23.92 

-18.58 

1 

B2 

0.86 

0.67 

-0.86 

-0.50 

1 

B3 

-38.27 

-36.95 

-6.39 

-2.23 

1 

B4 

-14.90 

-14.97 

-15.85 

-13.56 

2 

B5 

-3.83 

-0.34 

-6.14 

-1.09 

2 

B6 

-23.47 

-23.36 

-21.35 

-18.62 

3 

B7 

-19.20 

-19.87 

-9.50 

-8.82 

3 

B8 

-24.55 

-23.15 

-18.52 

-17.14 

4 

B9 

-19.29 

-21.31 

-18.84 

-14.66 

4 
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Table  6-6.  Summary  of  activity  estimation  errors  for  the  heterogeneous  density 
experiments  by  using  multiple  energy  peaks  in  the  CGNN  and  MLEM  algorithms 


Error  <50%        <30%         <25%         <20%  <10% 


Type  A  Heterogeneity  :  total  9  cases 

2-Energv  Reconstruction 

CGNN  9'  7  7  5  3 

(100%)*'    (77.8%)       (77.8%)       (55.6%)        (33.3%) 


MLEM  9  8  7 


J 


(100%)  (88.9%)  (77.8%)  (66.7%)  (33.3%) 

3 -Energy  Reconstruction 

CGNN            9                 8                 8                 6  2 

(100%)  (88.9%)  (88.9%)  (66.7%)  (22.2%) 

MLEM           9                  8                 8                8  3 

(100%)  (88.9%)  (88.9%)  (88.9%)  (33.3%) 

Error  <50%  <30%  <25%  <20%  <10% 


Type  B  Heterogeneity  :  total  9  cases 

2-Energv  Reconstruction 

CGNN             9               8                8  5                 2 

(100%)      (88.9%)      (88.9%)  (55.6%)       (22.2%) 

MLEM            9                9                9  7                 4 

(100%)      (100%)       (100%)  (77.8%)       (44.4%) 

3-Energv  Reconstruction 

CGNN              9                 8                 8  4                  2 

(100%)       (88.9%)      (88.9%)  (44.4%)       (22.2%) 

MLEM             9                 9                  9  9                  4 

(100%)        (100%)       (100%)  (100%)        (44.4%) 

Number  of  cases  within  the  specified  error  bound. 
Percentage  of  the  9  cases  within  the  specified  error  bound. 
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nonuniform  distribution  of  an  Ir'^^  source  in  a  computational  grid  caused  observed 
detector  responses  to  deviate  from  calculated  ideal  detector  responses.  Such  local 
nonuniformity  of  source  distribution  was  especially  important  to  the  modeling  of 
numerical  simulations  described  in  the  next  chapter.  The  following  paragraphs  discuss 
these  two  findings. 

Results  of  the  36  experiments  described  in  the  last  two  sections  reveal  that  the 
computational  model  using  multiple  energy  peaks  in  the  CGNN  or  MLEM  algorithms 
was  able  to  predict  the  total  activity  of  randomly  distributed  sources  inside  a  waste  barrel 
at  a  high  accuracy  level.  Table  6-7  summarizes  all  results.  All  experiments  were  within 
50%  accuracy.  Using  2-energy  reconstruction,  the  CGNN  predicted  32  cases  (88.9%) 
and  the  MLEM  predicted  34  cases  (94.4%)  with  30%  accuracy.  Using  3-energy 
reconstruction,  the  CGNN  predicted  34  cases  (94.4%)  and  the  MLEM  predicted  35  cases 
(97.2%).  With  an  accuracy  level  as  high  as  20%,  the  MLEM  predicted  30  cases  (83.3%) 
by  using  3-energy  reconstruction.  These  experimental  facts  verified  the  prediction  ability 
of  the  computer  model  developed  in  this  research.  However,  the  number  of  experiments 
was  still  too  limited  to  be  able  to  define  a  statistically  reasonable  confidence  level  for  the 
computer  model.  Further  and  more  extensive  verifications  were  done  by  using  the 
numerical  experiments  described  in  the  next  chapter. 

An  important  fact  needs  to  be  carefully  examined  in  regard  to  the  locally 
nonuniform  distribution  of  Ir'^^  sources  in  a  computational  grid  for  these  36  experiments. 
The  calculation  model  used  37  computational  grids  for  a  vertical  segment  of  the  waste 
barrel  and  assumed  that  once  a  source  was  in  a  grid,  the  source  was  uniformly  distributed 
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Table  6-7.  Summary  of  activity  prediction  for  36  experiments 


Error  <50%         <30%  <25%  <20%  <10% 


2-Energv  Reconstruction 

CGNN            36*  32  28  21  13 

(100%)'*  (88.9%)  (77.8%)  (58.3%)  (36.1%) 

MLEM          36  34  29  25  15 

(100%)  (94.4%)  (80.6%)  (69.4%)  (41.7%) 

3 -Energy  Reconstruction 

CGNN          36  34  29  22  12 

(100%)  (94.4%)  (80.6%)  (61.6%)  (33.3%) 

MLEM         36  35  32  30  17 

(100%)  (97.2%)  (88.9%)  (83.3%)  (47.2%) 


Number  of  cases  within  the  specified  error  bound. 
Percentage  of  the  36  cases  within  the  specified  error  bound. 
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within  that  grid.  It  is  difficult,  however,  to  simulate  a  fully  uniform  source  distribution 
condition  within  a  computational  grid  by  using  a  Ir'"-  source  of  only  3  mm  in  length  and 
0.5  mm  in  diameter.  Also,  in  the  experiments,  no  specific  bias  and  no  attempt  were  made 
to  place  Ir"^  pieces  in  acrylic-glass  compartments  such  that  a  uniform  distribution 
condition  could  be  simulated.  Therefore,  the  experimental  source  distribution  for  a 
computational  grid  deviated  somewhat  from  the  uniform  assumption.  Because  of  this 
local  nonuniformity,  observed  detector  responses  must  have  deviations  from  the  ideal 
calculation.  This  deviation  is  modeled  as  a  nonuniformity  noise  in  this  research.  In  the 
real  world,  radioactive  sources  inside  the  waste  drum  are  usually  clustered  and  rarely 
have  a  point  source  inside  a  computational  grid.  Local  nonuniformity  noise  is  more 
severe  in  a  point  source  case  than  in  a  clustered  source  case.  The  reason  for  this  is  that 
for  a  clustered  source,  observed  counts  of  a  detector  result  from  photon  emissions 
averaged  over  a  certain  volume  which  is  often  larger  than  the  size  of  the  Ir"^  source  used 
in  the  above  experiments.  Source  emissions  averaged  over  a  larger  volume  in  a 
computational  grid  have  a  greater  tendency  toward  a  uniform  emission,  meaning  that  the 
above  experiments  had  the  inherent  noise  caused  by  the  local  nonuniformity  of  the  small 
piece  of  Ir'^^  inside  a  grid.  The  magnitude  of  this  noise  may  compete  with  that  of  other 
noise  caused  by  statistics  of  the  detector  reading  or  the  uncertainties  of  the  calibration 
parameters.  Differences  of  magnitude  between  these  noises  are  analyzed  in  the  following 
paragraphs. 

First,  the  efficiency  calibration  was  very  trivial  and  is  easy  to  control  its 
uncertainty  below  a  few  percent.  For  example,  the  efficiency  of  3 1 6  keV  energy  peaks  of 
Ir^'^  were  found  to  be  2.17921x10^  and  this  energy  was  very  close  to  one  of  the 
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calibration  points  of  the  302.8  keV  of  Ba'''  check  source  (see  Table  6-1  and  Figure  6-3). 
The  efficiency  of  the  302.8  keV  gamma-ray  was  measured  as  2.2530x10-^  and  the  fitted 
efficiency  based  on  the  equation  of  the  calibration  curve  for  this  energy  was  2.2467x1 0"\ 
The  relative  deviation  was  only  0.28%.  The  uncertainty  (i.e.  the  noise)  of  efficiency 
calibration  for  the  3 16  keV  energy  peak  was  expected  to  be  close  to  this  value. 
Therefore,  the  contribution  of  noise  of  the  ideal  detector  response  from  the  efficiency 
calibration  was  very  small. 

Second,  the  magnitude  of  local  nonuniformity  noise  can  be  examined  by 
comparing  minimum,  maximum,  and  average  values  of  calculated  detector  kernel 
responses.  Detector  kernel  responses  for  the  homogeneous  condition  (listed  in  Appendix 
D)  were  calculated  by  the  computer  model  according  to  Equation  3- 1 .  The  3 1 6  keV 
kernel  responses  for  grid  19  or  grid  36  with  respect  to  detector  position  1  were 

min^,  ,9  =  2.9166x1 0""         maxK,,,  =  1.0995x10"'         oveK,,,  -  6.0387x10"' ; 

min/C,3g  =  1.9622x1 0-"         maxiC,  3,  =  1.2040x10"'         aveK,,,  =  5.3356x10"'; 
where  aveKjj  was  the  average  kernel  response  for  detector  position  i  contributed  by 
5,000  Monte  Carlo  samplings  of  point  sources  in  grid/;  and  minAT/y  was  the  minimum 
and  maxKij  was  the  maximum  point  kernel  response  among  the  5,000  points, 
respectively.  Deviations  from  the  average  for  these  two  grids  were 

'""""^'■"-""^'"^^  ^63.22%         """'*^'-"-°''^^'"L  12565% 


aveK^^g  aveK, 


1,19 
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A  similar  induction  used  for  grid  36  with  respect  to  detector  position  4  was 

(min^,3, -Qvg^  3^ )  (max  K, ,,  -  aveK, ,, ) 

— =  -57.54% ;       '-^ 'J^  =  1 05.1 0% 

These  calculations  reveal  that  from  the  point  of  view  of  detector  position  1,  a  point  source 
at  the  very  far  comer  of  grid  36  contributed  to  a  detector  response  which  was  63.22% 
lower  than  that  contributed  by  a  uniformly  distributed  source  in  the  grid,  and  a  point 
source  at  the  very  near  comer  contributed  to  a  detector  response  which  is  125.65%)  higher 
than  that  contributed  by  the  uniformly  distributed  source.  These  two  values  defined  the 
low  and  high  limits  of  the  local  nonuniformity  noise  of  grid  36  with  respect  to  detector 
position  1.  For  grid  19,  the  low  and  high  limits  of  noise  became  -51.70%)  and  82.08%. 
For  grid  36  with  respect  to  detector  position  4,  the  low  and  high  limits  of  noise  were 
-57.54%  and  105.10%.  From  these  comparisons,  it  is  evident  that  local  nonuniformity  of 
a  point  source  caused  a  wide  range  of  perturbation  of  the  ideal  detector  response  which 
was  calculated  based  on  the  average  kernel  response  of  each  computational  grid.  This 
situation  was  more  severe  in  a  peripheral  grid  than  in  a  central  grid.  Table  6-8  shows  the 
comparisons  of  ideal  and  observed  counting  rates  for  experiment  case  HI  (1  Ir'^^  piece  in 
central  grid  19)  and  case  H3  (1  Ir"'  piece  in  peripheral  grid  36).  These  two  cases 
represent  two  extremes  of  source  distribution  inside  the  waste  barrel.  Obviously,  case  H3 
was  much  more  perturbed  than  case  HI ,  as  shown  in  Table  6-8.  Case  HI  was  a  fully 
synmietric  case  and  all  four  detector  positions  were  ideally  supposed  to  have  the  same 
readings.  Case  H3  was  a  symmetric  case  for  detector  2  and  detector  4,  and  ideally  these 
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Table  6-8.  Comparisons  of  ideal  and  observed  detector  responses  for  experiment  case  HI 
and  case  H3 


Case  HI 

316  keV 

468  keV 

588  keV 

Ideal 

Obs. 

%Dev.* 

Ideal 

Obs. 

%Dev. 

Ideal    Obs. 

yoDev. 

Detr* 

13.60 

15.09 

10.96 

9.41 

10.31 

9.56 

0.97     1.06 

9.28 

Det2 

13.60 

13.95 

2.57 

9.41 

10.10 

7.33 

0.97     0.99 

2.06 

Det3 

13.60 

14.22 

4.56 

9.41 

9.84 

4.57 

0.97     0.93 

-4.12 

Det4 

13.60 

13.39 

-1.54 

9.41 

9.47 

0.64 

0.97     0.96 

-1.03 

Case  H3 

316  keV 

468  keV 

588  keV 

Ideal 

Obs. 

%Dev. 

Ideal 

Obs. 

%Dev. 

Ideal    Obs. 

%Dev. 

Detl   1190.00 

1024.00 

-13.95 

544.60 

481.10 

-11.66 

43.93    39.01 

-11.20 

Det2 

6.48 

7.19 

10.96 

4.71 

5.21 

10.62 

0.50     0.65 

30.00 

Det3 

0.34 

0.30 

-11.76 

0.35 

0.32 

-8.57 

0.04      0.02 

-50.00 

Det4 

6.48 

9.78 

50.93 

4.71 

6.76 

43.52 

0.50      0.68 

36.00 

Percent  of  deviation  of  observed  response  from  ideal  responses 
Detector  position  1 
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two  detectors  were  supposed  to  have  the  same  readings.  However,  no  obvious  symmetry 
was  observed  from  either  of  these  two  experiments. 

The  standard  deviation  of  the  counting  rate  for  the  counting  statistics  was 
(7  =  ^D/t  {D  being  counting  rate  and  t  being  total  counting  time).  For  a  300  second 
counting  time,  the  standard  deviation  of  the  3 16  kev  counting  rate  for  detector  position  1 
was  0.22  in  experiment  HI  and  was  1.85  in  experiment  H2.  These  were  only  1.62%  and 
0.16%  of  the  ideal  responses  for  these  two  cases  respectively  -  far  below  the  percentage 
deviations  of  10.96%  and  -13.95%  as  shown  in  Table  6-8.  Therefore,  counting  statistics 
was  not  the  major  factor  causing  the  observed  response  deviated  from  the  ideal  response 
for  these  two  experiments.  Perturbations  caused  by  local  nonuniform  source  distributions 
were  the  important  contributors  to  the  detector  response  noise  in  the  experiments. 


CHAPTER  7 
NUMERICAL  SIMULATIONS 

In  addition  to  36  experiments  described  in  Chapter  6,  this  research  also  performed 
numerical  simulations  to  determine  the  error  bound  and  its  associated  confidence  level  for 
the  computer  model  in  a  reasonable  statistical  sense.  In  real  experimentation,  the 
computer  model  uses  backward  reconstruction  algorithms  CGNN  and/or  MLEM  to 
estimate  total  activity  based  on  observed  external  detector  responses.  In  numerical 
experimentation,  total  activity  was  obtained  from  backward  reconstruction  based  on  ideal 
detector  responses  which  were  calculated  by  a  forward  projection  algorithm.  However, 
as  discussed  in  the  last  section  of  Chapter  6,  observed  responses  are  always  different  from 
ideal  responses.  These  differences  were  treated  as  "noises"  in  this  research.  To  simulate 
real  measurements  by  using  numerical  experiments,  "noise"  needed  to  be  modeled  and 
incorporated  in  ideal  responses.  This  research  modeled  two  types  of  noise  -  local 
nonuniformity  noise  and  counting  statistical  noise.  Since  waste  barrel  contents  are 
mostly  unknown,  density  heterogeneity  has  a  certain  degree  of  influence  on  the  accuracy 
of  computed  result.  This  research  used  a  random  density  perturbation  model  to  simulate 
different  types  of  density  heterogeneity  and  combined  these  density  perturbations  with 
noise  perturbations  to  form  a  wide  spectrum  of  numerical  simulations. 
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The  Density  Perturbation  Model 

Density  heterogeneity  was  modeled  by  specifying  a  relative  density  perturbation 
level  with  respect  to  a  homogeneous  base  density.  The  perturbed  density  for  each 
computational  grid  was  assigned  by  the  Monte  Carlo  process: 

P;=p"-[l  +  /-(2/7y-l)]  (7-1) 

where  p'^  was  the  perturbed  density  for  the  computational  grid;,  p°  was  the  base 
density,  /  was  the  relative  density  perturbation  level,  and  7y  was  the  random  number 
generated  for  the  computational  grid  7. 

In  order  to  have  a  consistent  basis  of  comparison  with  experiments  described  in 
Chapter  6,  1.0  g/cm^  was  used  as  the  base  density  in  the  simulation.  Two  relative  density 
perturbation  levels,  20%  (low)  and  50%  (high),  were  studied.  UncoUided  gamma-rays 
emitted  from  any  randomly  sampled  source  point  needed  to  pass  one  or  several  grids  to 
reach  the  detector.  Path  length  within  each  grid  was  accurately  calculated  for  density 
perturbation  conditions.  For  each  simulated  case,  a  perturbed  kernel  table  was  obtained 
by  using  Equation  3-1,  based  on  calculated  path  length  and  density  distribution  calculated 
from  Equation  7-1 .  Then,  the  kernel  table  was  incorporated  with  the  assignment  of 
random  source  fraction  of  each  grid,  to  give  the  forward  projection  of  ideal  detector 
responses  by  using  Equation  3-2. 

The  Model  of  Local  Nonuniformity  Noise 

Ideal  detector  responses  were  generated  using  the  assumption  that  the  source  was 
uniformly  distributed  within  a  computational  grid.  As  found  in  Chapter  6,  local 
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nonuniformity  of  source  distribution  within  a  computational  grid  was  an  important 
contributor  of  signal  noise  and  was  modeled  by  specifying  a  relative  perturbation  level 
for  detector  responses.  This  level  was  the  perturbation  range  of  ideal  detector  responses, 
which  was  assumed  to  occur  in  real  experiments.  A  Monte  Carlo  process  was  used  to 
select  the  nonuniformity  noise  from  this  range.  That  is,  for  each  detector  position  i,  the 
perturbed  response  (counting  rate)  was 

D;iE^)  =  Df(E^)-[l  +  p.(2-rj,-\)]  (7-2) 

where  D;(Ej  )  was  the  perturbed  response  for  energy  peak  j  of  detector  position  /; 
D'{Ej)  was  the  ideal  detector  response;  p  was  the  relative  perturbation  level  contributed 
from  local  nonuniformity  noise;  and  r/.  was  the  random  number  generated  for  detector 
position  /.  Perturbation  levels  from  0%  to  50%  were  studied.  As  shown  in  Table  6-8,  the 
50%  perturbation  case  was  equivalent  to  experiment  case  H3  which  was  a  highly  biased 
case  with  an  Ir''-  source  very  close  to  one  of  the  detector  positions  and  having  a 
maximum  deviation  around  50%  between  observed  and  ideal  responses. 

The  Model  of  Counting  Statistics 

Any  radiation  measurement  is  subject  to  the  uncertainty  of  counting  statistics.  A 
perturbed  detector  response  £),'(£-)  calculated  from  Equation  7-2  was  considered  an 
average  value  of  the  measurement  of  randomly  distributed  sources  nonuniformly 
distributed  inside  any  computational  grid.  To  simulate  a  real  response  with  counting 
statistical  behavior,  distribution  of  perturbed  detector  response  D/(£^ )  was  modeled  by 
using  a  Gaussian  distribution  form 
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PiX)  =  -j^.,^^C^^-p')  (7.3) 


I  r\'(  p   \ 

where  X=  D;{Ej),  cr  =  J    '      '     ,  and  7  was  the  counting  time.  Using  the  rejection 

technique,  each  detector  response  in  the  numerical  simulation  was  obtained  by  sampling 
from  Equation  7-3,  with  a  sampling  range  from  J  -  Bcr  to  X  +  3cr .-' 

Error  Bound  and  Confidence  Level 

Five  numerical  experiment  categories  were  designed  with  the  following 
simulation  conditions: 

(i)   homogeneous  and  noise-free, 

(ii)  heterogeneous  with  20%  density  perturbation  and  noise-free, 
(iii)  heterogeneous  with  50%  density  perturbation  and  noise-free, 
(iv)  heterogeneous  with  20%  density  perturbation  and  different  noise  levels,  and 
(v)    heterogeneous  with  50%  density  perturbation  and  different  noise  levels. 
For  each  category,  evaluations  of  error  bound  and  its  associated  confidence  level 
for  activity  estimation  were  carried  out  by  using  different  backward  reconstruction 
schemes  under  different  simulation  conditions.  Each  simulation  is  performed  with 
respect  to  10,000  random  source  distribution  cases.  Confidence  level  for  the  simulation 
was  defined  as 
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with  Ai  being  the  estimated  activity  for  any  single  case,  A^  the  reference  activity,  e  the 
required  accuracy,  and  A^(=  10,000)  the  total  number  of  cases.  P^  is  the  frequency  of 
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occurrence  for  relative  error  within  £• .  For  example,  if  the  error  bounds  is  specified  as 
10%,  and  95%  of  the  predicted  cases  fell  within  10%  of  the  actual  source  strength,  the 
prediction  error  would  be  said  to  be  10%  at  the  95%  confidence  level. 

In  the  simulation,  counting  time  was  300  seconds  and  reference  activity  was 
10^  Bq.  This  activity  level  was  about  the  same  order  of  magnitude  as  those  activities  used 
in  the  experiments  described  in  Chapter  6. 

The  Naming  System  for  the  Numerical  Experiment 

In  order  to  distinguish  the  results  of  a  large  number  of  simulation  cases,  a  naming 
system  was  designed.  For  category  (i),  six  evaluations  were  studied  and  were  named  as 
HCEl,  HCE2,  HCE3,  HMEl,  HME2,  and  HME3.  "H"  stands  for  homogeneous,  "C" 
stands  for  the  CGNN  method,  "M"  stands  for  the  MLEM  method,  and  "El",  "E2",  and 
"E3"  denote  1 -energy,  2-energy  and  3 -energy  reconstruction,  respectively.  For  categories 
(ii)  to  (v),  the  naming  system  used  the  notation  "XnNmmYEp"  explained  by  the 
following  keys: 

"X"  is  "C"  or  "M"  and  represents  the  CGNN  or  the  MLEM  method,  respectively; 

"n"  is  "2"  or  "5"  and  means  20%  or  50%  density  perturbation,  respectively; 

"N"  denotes  local  nonuniformity  noise; 

"mm"  is  the  relative  perturbation  level  of  local  nonuniformity  noise; 

"Y"  is  "B"  or  "K"  and  represent  the  reconstruction  is  based  on  the  base  density 
kernel  table  or  the  reconstruction  is  based  on  the  perUirbed  kernel  table,  respectively; 

"E"  means  energy;  and 
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"p"  is  "1",  "2"  or  "3"  and  represents  1-energy,  2-energy,  and  3-energy 
reconstruction,  respectively. 

For  example,  a  simulation  result  denoted  by  C2N05KE2  would  mean  that  the 
result  was  based  on  using  2  energies  and  the  perturbed  kernel  table  in  the  CGNN 
reconstruction  algorithm  under  a  20%  density  perturbation  and  a  5%  noise  condition.  A 
result  denoted  by  M5N50BE2  would  mean  that  the  result  was  based  on  using  2  energies 
and  the  base  density  kernel  table  in  the  MLEM  reconstruction  algorithm  under  a  50% 
density  perturbation  and  a  50%)  noise  condition. 

Homogeneous  and  Noise-Free  Condition 

Table  7-1  shows  simulation  results  of  different  reconstruction  schemes  under  the 
homogeneous  and  noise-free  condition  (simulation  category  (i)).  In  this  simulation 
category,  1-energy,  2-energy,  and  3-energy  reconstructions  were  used  for  both  the  CGNN 
and  MLEM  algorithms.  Confidence  levels  of  1-energy  reconstruction  were  very  low. 
However,  when  two  energies  were  used,  confidence  levels  dramatically  increased  to  a 
very  satisfactory  result  for  both  CGNN  and  MLEM  reconstructions.  All  the  predictions 
of  the  CGNN  were  within  20%)  error  and  the  confidence  level  reached  98.1%)  for  an 
accuracy  as  high  as  10%.  All  the  predictions  of  the  MLEM  were  within  25%  error  and 
the  confidence  level  for  10%o  error  was  93.5%).  Confidence  level  increased  even  more 
when  three  energies  were  used,  but  not  as  dramatically  as  when  increasing  from  one 
energy  to  two  energies.  Almost  all  the  predictions  (99.7%))  were  within  10%o  error  for 
3-energy  CGNN  reconstruction.  The  confidence  level  for  10%  error  was  92.1%  for 
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Table  7-1.  Confidence  levels  for  total  activity  estimation  by  different  optimization 
algorithms  under  the  homogeneous  and  noise-free  condition 


Category  (i) :  Homogeneous  and  Noise-Free  Condition 


Error  <50%  <30%  <25%        <20%  <10% 

HCEl*  0.438  0.301  0.267  0.229  0.120 

HMEl"  0.630  0.372  0.323  0.277  0.158 

HCE2  1.000  1.000  1.000  1.000  0.981 

HME2  1.000  1.000  1.000  0.994  0.935 

HCE3  1.000  1.000  1.000  1.000  0.997 

HME3  1.000  1.000  1.000  1.000  0.921 


*^  H:  homogeneous  density;  C:  CGNN;  M:  MLEM; 
"El:  1-energy;  E2:  2-energy;  E3:  3 -energy. 
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3-energy  MLEM  reconstruction.  In  this  simulation  condition  (category  (i)),  prediction 
ability  of  the  CGNN  proved  better  than  that  of  the  MLEM. 

Figure  7-1  shows  the  distributions  of  relative  prediction  error  for  2-energy  CGNN 
and  MLEM  methods  in  simulation  category  (i).  It  is  even  more  evident  that  the  CGNN's 
predictions  were  all  within  15%  error  and  had  no  significant  bias  toward  underestimation 
or  overestimation.  For  the  MLEM,  a  tendency  toward  overestimation  can  be  seen,  but  the 
error  is  not  more  than  25%. 

Comparisons  with  the  Segmented  Gamma  Scanning  Method 

The  traditional  Segmented  Gamma  Scanning  (SGS)  method  assumes  that  the 
material  and  radionuclides  are  distributed  homogeneously  and  uniformly  within  the  waste 
drum.  Data  collection  is  done  by  measuring  the  gamma  signal  emitted  from  the  drum  as 
it  is  rotated  stepwisely  or  continuously.  The  rotation  is  used  to  compensate  for  activity 
distribution  asymmetries.  Under  this  assumption,  the  total  detector  response  DT  is  the 
sum  of  all  detector  responses  and  is  proportional  to  total  activity: 

i 

For  1  Bq  homogeneous  source,  the  total  response  DThom  can  be  easily  calculated  by 

doing  an  analytical  integration  or  a  Monte  Carlo  integration  over  the  drum  volume. 

Therefore,  for  any  measured  detector  response  Dr^^n  which  is  essentially 

nonhomogeneous  in  nature,  total  activity  is  estimated  to  be 

DT 
5,.  =— f^     (Bq). 


67 


u 
o 

O 


cr 


-65   -55   -45   -35   -25   -15-5      5      15     25     35     45     55     65 

-70-60-50^-30-20-10  0   10  20  30  40  50  60  70 
Relative  Error  of  Total  Activity  Estimation  (%) 


o 

0) 


o 

15 

cr 

0) 


0  50 

J 1 \ 1 \ 1 L 

L 

_! 

___.!      1      1      1      1 

0.45     ~ 

- 

0.40     - 

3 

Case  -  HME2 

= 

\ 

~ 

0.35     - 

^ 

= 

\ 

— 

\ 



0.30     - 

\ 

— 

0.25     - 

^ 

~ 

N 

\ 

_ 

\ 

\ 

— 

0.20     - 

^ 

^ 

— 

0.15     - 

^ 

1 

— 

0.10     ~ 

^ 

^   N 

^ 

X 

— 

0.05     - 
000 

J—  _i 

^ 

^ 

^ 

- 

n    1    \    \    \    1    1 

r^T    \    \ — \ — r 

-65    -55    -45    -35    -25    -15-5       5      15     25     35     45     55     65 

-70-60-50^-30-20-10  0   10  20  30  40  50  60  70 
Relative  Error  of  Total  Activity  Estimation  (%) 


Figure  7-1 .  Distributions  of  relative  prediction  error  for  the  2-energy  CGNN  and  MLEM 
algorithms  under  the  homogeneous  and  noise-free  condition 
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Table  7-2  shows  the  results  of  the  SGS  method  using  various  multiple  detector 
positions  and  compares  them  with  the  results  of  the  2-energy  CGNN  and  MLEM 
methods.  It  is  evident  that  increasing  the  number  of  detector  positions  in  SGS  did  not 
improve  the  accuracy  of  the  measurement.  For  prediction  within  30%  error,  the 
confidence  level  was  34.3%  at  best.  Even  for  an  error  bound  as  large  as  50%,  the 
confidence  level  in  the  SGS  measurement  was  very  low  and  less  than  d&'A.  Distributions 
of  relative  prediction  error  for  4  and  64  detector  positions  are  shown  in  Figure  7-2.  The 
maximum  error  for  4  detector  positions  occurred  at  +170%.  For  64  detector  positions, 
the  maximum  error  occurred  at  +100%.  Compared  with  the  2-energy  CGNN  or  MLEM 
method,  the  SGS's  prediction  was  very  poor. 

Heterogeneous  and  Noise-Free  Conditions 

Simulation  results  of  2-energy  and  3-energy  CGNN  and  MLEM  algorithms  under 
heterogeneous  and  noise-free  conditions  (simulation  categories  (ii)  and  (iii) )  are  shown  in 
Table  7-3.  For  20%)  heterogeneity,  2-energy  and  3-energy  CGNN  reconstruction  schemes 
predicted  10%  error  at  97%  and  98%)  confidence  levels  respectively.  Under  the  same 
conditions  for  the  MLEM,  2-energy  and  3-energy  reconstructions  predicted  10%)  error  at 
91.2%  and  91.5%  confidence  levels  respectively.  For  50%  heterogeneity,  the  CGNN 
predicted  10%)  error  with  92.9%o  confidence  by  using  2  energies  and  10%  error  with 
96.2%  confidence  by  using  3  energies.  The  MLEM  predicted  20%o  error  at  a  confidence 
level  of  less  than  90%,  but,  for  25%  error,  the  predictions  had  a  confidence  level  of  more 
than  92%.  In  these  two  categories,  3-energy  reconstructions  performed  better  than 
2-energy  reconstructions,  and  the  CGNN's  prediction  was  better  than  that  of  MLEM. 


69 


Table  7-2.  Comparisons  of  confidence  levels  between  the  SGS  and  the  2-energy  CGNN 
and  MLEM  reconstruction  algorithms  under  the  homogeneous  and  noise-free  condition 


Error  <50% 

<30% 

<25% 

<20% 

<10% 

SGS04* 

0.589 

0.277 

0.226 

0.174 

0.093 

SGS08' 

0.653 

0.343 

0.266 

0.212 

0.100 

SGS  16' 

0.656 

0.336 

0.278 

0.218 

0.105 

SGS32' 

0.657 

0.335 

0.278 

0.217 

0.105 

SGS64* 

0.657 

0.336 

0.279 

0.217 

0.105 

HCE2*' 

1.000 

1.000 

1.000 

1.000 

0.981 

HME2*" 

1.000 

1.000 

1.000 

0.994 

0.935 

Digits  indicate  the  number  of  detector  positions  for  SGS  method. 
Homogeneous  density  and  2-energy  CGNN  reconstruction. 
Homogeneous  density  and  2-energy  MLEM  reconstruction. 
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Figure  7-2.  Distributions  of  relative  prediction  error  for  the  SGS  method  using  4  and  64 
detector  positions  under  the  homogeneous  and  noise-free  condition 
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Table  7-3.  Confidence  levels  for  total  activity  estimation  by  different  optimization 
algorithms  under  the  heterogeneous  and  noise-free  condition 

Category  (ii) :  20%  Heterogeneity  and  Noise-Free  Condition 

Error  <50%        <30%         <25%         <20%        <10% 

C2N00KE2*  1.000  1.000  1.000  1.000  0.970 

M2N00KE2  1.000  1.000  0.999  0.995  0.912 


C2N00KE3 

1.000 

1.000 

1.000 

1.000 

0.980 

M2N00KE3 

1.000 

1.000 

1.000 

0.998 

0.915 

Category  (iii)  :  50%  Heterogeneity  and  Noise-Free  Condition 


Error   <50%       <30%        <25%         <20%        <10% 


C5N00KE2 

1.000 

1.000 

0.999 

0.996 

0.929 

M5N00KE2 

0.997 

0.953 

0.923 

0.851 

0.599 

C5N00KE3 

1.000 

1.000 

0.998 

0.998 

0.962 

M5N00KE3 

0.999 

0.967 

0.939 

0.886 

0.616 

*  Notation  for  the  1st  column-  C2:  CGNN  and  20%  heterogeneity; 
M2:  MLEM  and  20%  heterogeneity;  C5:  CGNN  and  50%  heterogeneity; 
M5:  MLEM  and  50%  heterogeneity;  NOO:  noise-free;  E2:  2-energy; 
E3:  3 -energy;  K:  perturbed  density  kernel  table. 
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Distribution  of  relative  prediction  error  for  2-energy  CGNN  and  MLEM 
algorithms  are  shown  in  Figure  7-3.  Maximum  relative  enror  for  the  2-energy  MLEM 
method  was  +30%  and  +60%  for  20%  and  50%  heterogeneity,  respectively.  For  the 
2-energy  CGNN  method,  maximum  relative  errors  was  -20%  and  -30%  for  20%  and  50% 
heterogeneity,  respectively. 

Heterogeneous  and  Noise-Perturbed  Conditions 

Local  nonuniformity  noise  and  counting  statistics  were  both  simulated  under 
heterogeneous  conditions.  Five  relative  perturbation  levels  of  local  nonuniformity  (1%, 
5%,  10%),  20%)  and  50%)  were  simulated.  Gaussian  statistics  were  applied  for  every 
noise-perturbed  detector  response. 
Simulation  Resuhs  of  the  20%  Heterogeneity  Level 

Table  7-4  shows  simulation  results  for  different  noise  levels  under  20%  density 
heterogeneity  (simulation  category  (iv) ).  The  results  of  1 -energy,  2-energy,  and  3-energy 
CGNN  and  MLEM  reconstruction  schemes  are  compared  at  5%  noise  level.  It  is  evident 
again  that  the  total  activity  predictions  of  1 -energy  CGNN  and  MLEM  methods  were 
very  poor.  A  surprising  finding  is  that  CGNN  results  deteriorated  by  changing  from  a  2- 
energy  to  a  3-energy  reconstruction  scheme.  Confidence  levels  for  30%  accuracy 
decreased  from  96.4%  to  93.7%,  and  for  25%  accuracy,  decreased  from  92.8%  to  89.7%. 
This  deterioration  phenomenon  was  not  found  in  the  3-energy  MLEM  results  at  this  noise 
level. 

Table  7-4  also  clearly  shows  influences  of  noise  at  this  heterogeneity  level.  For 
2-energy  CGNN  reconstruction,  the  confidence  level  for  30%)  accuracy  decreased  from 
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Figure  7-3.  Distributions  of  relative  prediction  error  for  2-energy  CGNN  and  MLEM 
methods  under  the  heterogeneous  and  noise-free  condition 
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Table  7-4.  Confidence  levels  for  total  activity  estimation  by  different  optimization 
algorithms  under  the  20%  heterogeneity  and  noise-perturbed  condition 

Category  (iv)  :  20%  Heterogeneity  with  Local 

Nonuniformitv  and  Statistical  Noise 


Error  <50%        <30%        <25%        <20%       <10% 


C2N01KE2 

0.998 

0.969 

0.932 

0.860 

0.563 

M2N01KE2 

0.998 

0.968 

0.944 

0.877 

0.600 

C2N05KE1 

0.468 

0.319 

0.278 

0.237 

0.141 

M2N05KE1 

0.663 

0.395 

0.337 

0.282 

0.160 

C2N05KE2 

0.998 

0.964 

0.928 

0.857 

0.557 

M2N05KE2 

0.998 

0.970 

0.944 

0.866 

0.588 

C2N05KE3 

0.993 

0.937 

0.897 

0.831 

0.553 

M2N05KE3 

0.999 

0.973 

0.951 

0.894 

0.638 

C2N10KE2 

0.999 

0.963 

0.919 

0.850 

0.529 

M2N10KE2 

0.999 

0.968 

0.940 

0.865 

0.573 

C2N20KE2 

0.998 

0.951 

0.896 

0.805 

0.495 

M2N20KE2 

0.999 

0.965 

0.912 

0.826 

0.524 

C2N50KE2 

0.979 

0.821 

0.743 

0.630 

0.314 

M2N50KE2 

0.990 

0.840 

0.757 

0.639 

0.328 

Notation  for  the  1st  column-  C2:  CGNN  and  20%  heterogeneity; 
M2:  MLEM  and  20%)  heterogeneity;  E2:  2-energy;  E3:  3 -energy 
K:  perturbed  density  kernel  table;  NOl :  \%  noise;  N05:  5%  noise; 
NIO:  10%  noise;  N20:  20%  noise;  N50:  50%  noise. 
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96.9%  to  82.1%  when  the  noise  level  increased  from  1%  to  50%.  For  2-energy  MLEM 
reconstruction,  the  confidence  level  for  30%  accuracy  decreased  from  96.8%  to  84.0% 
when  the  noise  level  increased  from  1%  to  50%.  Compared  with  the  noise-free  cases 
described  in  the  last  two  sections,  the  performance  of  the  CGNN  is  reversed  with  respect 
to  the  MLEM  in  the  noisy  condition.  The  CGNN  confidence  level  is  degraded  faster  than 
that  of  the  MLEM  when  noise  increased.  Even  in  only  a  1%  noise  condition,  the  MLEM 
outperformed  the  CGNN.  This  phenomenon  can  be  explained  by  the  fact  that  the  MLEM 
is  theoretically  based  on  a  statistical  derivation  and  the  noise  is  also  modeled  statistically 
by  this  research.  Therefore,  the  MLEM's  performance  is  somehow  more  regulative  in  the 
noise-perturbed  condition. 

For  50%>  accuracy,  the  confidence  level  was  97.9%.  for  the  2-energy  CGNN 
method  and  99.0%.  for  the  2-energy  MLEM  method.  A  most  likely  accuracy  for  both  the 
CGNN  and  MLEM  methods  under  20%  heterogeneity  condition  can  be  determined  from 
Table  7-4.  It  is  defined  as  the  highest  accuracy  with  more  than  95%.  confidence  level 
which  can  be  found  from  the  table.  For  both  2-energy  CGNN  and  MLEM  results  shown 
in  Table  7-4,  no  confidence  level  greater  than  95%.  could  be  found  for  accuracy  higher 
than  30%.  The  confidence  level  for  30%  accuracy  was  95.1%  for  the  2-energy  CGNN 
method  and  96.5%  for  the  2-energy  MLEM  method  at  20%.  noise  level.  Therefore,  the 
most  likely  accuracy  for  both  of  them  was  said  to  be  30%  error  with  a  noise  level  up  to 
20%. 

Distributions  of  relative  prediction  error  for  2-energy  CGNN  and  MLEM  methods 
under  5%,  20%  and  50%.  noise  level  are  shown  in  Figure  7-4.  Maximum  relative  error 
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Figure  7-4.  Distributions  of  relative  prediction  eixor  for  2-energy  CGNN  and  MLEM 
methods  under  the  20%  heterogeneity  and  different  noise-perturbed  condition 
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for  the  CGNN  in  these  three  noise  levels  was  +60%,  +60%,  and  +70%,  respectively.  For 
the  MLEM,  it  was  +55%,  +60%,  and  +65%,  respectively.  The  frequency  of  occurrence 
for  all  these  maximum  relative  errors  was  below  0.5%). 
Simulation  Resuhs  of  the  50%  Heteroeeneitv  Level 

Table  7-5  shows  simulation  results  for  different  noise  levels  under  50%  density 
heterogeneity  (simulation  category  (v) ).  Figure  7-5  shows  the  distributions  of  relative 
prediction  error  for  2-energy  CGNN  and  MLEM  methods  under  5%,  20%  and  50%  noise 
levels.  The  results  of  2-energy  and  3 -energy  CGNN  and  MLEM  reconstruction  schemes 
are  compared  at  the  50%  noise  level  in  Table  7-5.  Again,  it  shows  that  for  the  CGNN, 
3 -energy  reconstruction  deteriorated  the  results  of  2-energy  reconstruction.  The 
confidence  level  for  50%  accuracy  decreased  from  97.8%  to  96.2%.  For  30%  accuracy, 
the  confidence  level  decreased  from  81.8%  to  78.6%.  For  the  MLEM,  the  prediction 
ability  of  3-energy  reconstruction  was  comparable  to  that  of  2-energy  reconstruction. 
Confidence  levels  for  30%  and  25%  accuracy  deteriorated,  but  for  20%>  and  10% 
accuracy,  the  confidence  level  is  improved.  This  means  that  the  overall  performance  of 
3-energy  MLEM  reconstruction  was  limited  when  the  noise  level  was  as  high  as  50%. 
Therefore,  for  both  CGNN  and  MLEM  application  in  the  high  noise  condition,  the 
2-energy  reconstruction  scheme  was  the  optimal  choice. 

For  50%  accuracy,  both  2-energy  CGNN  and  2-energy  MLEM  reconstruction 
schemes  provided  more  than  a  95%  confidence  level  under  a  noise  level  as  high  as  50%. 
These  results  verified  the  fact  that  36  experimental  results  described  in  Chapter  6  were  all 
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Table  7-5.  Confidence  levels  for  total  activity  estimation  by  different  optimization 
algorithms  under  the  50%  heterogeneity  and  noise-perturbed  condition 


Category  (v)  :  50%  Heterogeneity  with  Local 

Nonuniformity  and  Statistical  Noise 

Error  <50% 

<30% 

<25% 

<20% 

<10% 

C5N01KE2'  0.996 
M5N01KE2  0.997 

0.953 
0.953 

0.922 
0.919 

0.851 
0.853 

0.550 
0.585 

C5N05KE2    0.998 
M5N05KE2  0.996 

0.954 
0.953 

0.918 
0.916 

0.844 
0.850 

0.531 
0.573 

C5N10KE2    0.998 
M5N10KE2  0.996 

0.955 
0.956 

0.910 
0.910 

0.834 
0.840 

0.515 
0.554 

C5N20KE2    0.998 
M5N20KE2  0.997 

0.942 
0.951 

0.892 
0.893 

0.789 
0.814 

0.458 
0.498 

C5N50KE2    0.978 
M5N50KE2  0.986 
C5N50KE3    0.962 
M5N50KE3  0.986 

0.818 
0.828 
0.786 
0.825 

0.719 
0.738 
0.706 
0.730 

0.605 
0.628 
0.582 
0.636 

0.315 
0.319 
0.296 
0.335 

Notation  for  the  1st  column-  C5:  CGNN  and  50%  heterogeneity; 
M5:MLEM  and  50%  heterogeneity;  E2:  2-energy;  E3:  3 -energy 
K:  perturbed  density  kernel  table;  NOl:  1%,  noise;  N05:  5%  noise; 
NIO:  10%  noise;  N20:  20%  noise;  N50:  50%  noise. 
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Figure  7-5.  Distributions  of  relative  prediction  error  for  2-energy  CGNN  and  MLEM 
methods  under  the  50%  heterogeneity  and  different  noise-perturbed  condition 
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within  50%  error.  Using  the  2-energy  CGNN  reconstruction  scheme,  the  most  hkely 
accuracy  was  30%  error  with  a  noise  level  up  to  10%.  However,  using  2-energy  MLEM 
reconstruction  scheme,  the  most  likely  accuracy  was  found  to  be  30%  error  with  a  noise 
level  up  to  20%. 

For  the  CGNN,  the  maximum  relative  error  of  the  three  noise  levels  shown  in 
Figure  7-5  was  +60%,  +70%,  and  +75%,  respectively.  For  the  MLEM,  it  was  +65%, 
+65%,  and  +70%,  respectively.  The  frequency  of  occurrence  for  all  these  maximum 
relative  errors  was  below  0.5%. 

Simulations  of  Biased  Source  Distribution  Cases 

Using  the  Monte  Carlo  process  in  the  computer  model,  a  source  fi-action  for  each 
computational  grid  was  randomly  assigned.  However,  this  process  was  not  be  able  to 
generate  biased  source  distribution  cases.  In  order  to  check  if  these  source  distribution 
cases  could  produce  biased  total  activity  estimation,  a  total  of  100  biased  source 
distribution  cases  were  simulated  and  the  total  activity  of  each  of  them  were 
reconstructed  by  using  2-energy  CGNN  method.  These  100  cases  encompassed  six  types 
of  biased  source  distribution:  (a)  single  source  concentrated  in  a  grid,  (b)  clustered  source 
concentrated  in  a  few  peripheral  grids,  (c)  clustered  source  concentrated  in  a  few  central 
grids,  (d)  sources  distributed  symmetrically  in  peripheral  grids,  (e)  sources  scattered  in  a 
quadrant  or  half  of  the  37  grids,  and  (f)  sources  distributed  symmetrically  in  central  grids. 
Figure  7-6  shows  representative  configurations  for  these  six  types  of  source  distribution. 
Table  7-6  shows  that  resulting  confidence  levels  based  on  these  100  cases  were 
comparable  to  the  results  of  the  2-energy  CGNN  method  which  were  based  on  10,000 


81 


y' 

■^ 

N 

\ 

r333 

— ^ 

1.333 

1333 

. 

7 

\ 

_^ 

7 

(a)  Single  Source  Concentrated  in  a  Grid  (b)  Clustered  Sources  in  Peripheral  Grids 


y" 

\ 

/ 

' 

,333 

.333 

.333 

\ 

Vs| 

y 

f-r 

K 

/ 

.250 

250|\ 

t 

r" 

r_ 

i 

\     .250 

.250 

7 

\| 

^ 

\}y 

(c)  Clustered  Sources  in  Central  Grids  (d)  Sources  Distributed  Symmetrically  in  Peripheral  Grids 


^ 

** 

.001 

.001 

\ 

/ 

.010 

.010 

.001 

\ 

t- 

.500 

.200 

.010 

—*• 

.001 

n 

.146 

.100 

.010 

.001 

[~ 

. 

\ 

7 

x 

^ 

^ 

/ 

(e)  Sources  Scattered  in  a  Quadrant       (i)  Sources  Distributed  Symmetrically  in  Central  Grids 


Figure  7-6.  Representative  configurations  for  six  types  of  biased  source  distribution 
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Table  7-6.  Comparison  of  the  confidence  levels  of  10,000  random  source  distribution 
cases  and  100  biased  source  distribution  case 


Error  <50% 

<30% 

<25% 

<20% 

<10% 

C2N05KE2* 

0.998 

0.964 

0.928 

0.857 

0.557 

C2N05KE2-BS 

1.000 

0.960 

0.920 

0.860 

0.560 

C2N20KE2 

0.998 

0.951 

0.896 

0.805 

0.495 

C2N20KE2-BS 

1.000 

0.940 

0.880 

0.810 

0.500 

C2N50KE2 

0.979 

0.821 

0.743 

0.630 

0.314 

C2N50KE2-BS 

0.970 

0.810 

0.730 

0.640 

0.320 

C5N05KE2 

0.998 

0.954 

0.918 

0.844 

0.531 

C5N05KE2-BS 

1.000 

0.960 

0.930 

0.870 

0.570 

C5N20KE2 

0.998 

0.942 

0.892 

0.789 

0.458 

C5N20KE2-BS 

1.000 

0.940 

0.910 

0.810 

0.410 

C5N50KE2 

0.978 

0.818 

0.719 

0.605 

0.315 

C5N50KE2-BS 

0.990 

0.850 

0.760 

0.610 

0.290 

*  Notation  for  the  1st  column-  C2:  CGNN  and  20%  heterogeneity; 
C5:  CGNN  and  50%  heterogeneity;  E2:  2-energy; 
K:  perturbed  density  kernel  table;  BS:  biased  source  distribution  cases; 
N05:  5%  noise;  N20:  20%  noise;  N50:  50%  noise. 
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random  source  distribution  cases.  Distributions  of  relative  prediction  error  are  shown  in 
Figure  7-7.  It  is  evident  that  these  100  highly  biased  source  distribution  cases  do  not 
deteriorate  the  results  or  cause  significant  biased  effects  with  respect  to  those  presented  in 
the  previous  section. 

Simulations  Using  the  Base  Density  Kernel  Table 

The  numerical  simulations  of  the  heterogeneous  condition  presented  above  were 
performed  by  simulating  a  perturbed  density  distribution  and  constructing  a  perturbed 
kernel  table.  Using  the  perturbed  kernel  table  and  a  random  source  distribution,  the 
forward  projection  of  ideal  detector  responses  was  obtained.  Ideal  detector  responses 
were  perturbed  by  the  noise  effect  to  simulate  observed  detector  responses.  Total  activity 
was  then  reconstructed  based  on  the  known  perturbed  kernel  table  and  the  simulated 
detector  responses.  This  simulates  a  situation  in  which  transmission  tomography  or 
radiography  is  performed  and  detailed  density  information  is  provided.  In  a  situation  in 
which  detailed  density  information  can  not  be  provided,  for  instance  when  application  of 
transmission  tomography  is  impossible,  backward  reconstruction  of  total  activity  can  only 
be  implemented  by  using  the  appropriate  base  density  kernel  table.  A  base  density  kernel 
table  is  constructed  based  on  the  average  density  of  the  waste  barrel.  A  prior  knowledge 
of  the  weight  and  volume  of  the  barrel  or  a  transmission  density  measurement  can  help  on 
the  determination  of  the  average  density.  However,  resulting  confidence  levels  will  be 
different  from  results  based  on  a  known  density  distribution.  In  this  research,  perturbed 
density  was  calculated  from  Equation  7-1  in  which  the  base  density  used  is  1.0  g/cm^ 
Naturally,  the  simulation  system  had  an  average  density  close  to  1.0  g/cm^  as  well.  This 
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Figure  7-7.  Distributions  of  relative  prediction  error  for  2-energy  CGNN  method  using 
100  biased  source  distribution  cases  under  heterogeneity  and  noise-perturbed  conditions 


85 

section  presents  numerical  simulations  for  the  20%  and  50%  heterogeneity  conditions  in 
which  detailed  density  information  was  unknown  and  therefore,  the  base  density  kernel 
table  was  used  for  backward  reconstruction. 

Table  7-7  shows  the  comparisons  of  confidence  levels  for  the  prediction  of  total 
activity  using  perturbed  kernel  and  base  density  kernel  tables  under  the  20% 
heterogeneity  condition.  Table  7-8  shows  the  same  comparisons  under  the  50% 
heterogeneity  condition.  From  these  two  tables,  it  is  evident  that  accuracy  was  degraded 
when  using  the  base  density  kernel  table  as  compared  to  when  using  the  perturbed  kernel 
table.  This  is  because  density  heterogeneity  became  an  extra  noise  component  added  to 
ideal  detector  responses  which  were  calculated  using  the  perturbed  kernel  table.  Under 
20%  heterogeneity  (Table  7-7),  the  confidence  levels  of  30%  accuracy  for  the  CGNN 
were  91.1%  for  the  5%  noise  level,  87.7%  for  the  20%  noise  level,  and  76.1%  for  the 
50%.  noise  level.  For  the  MLEM,  these  confidence  levels  were  93.9%,  90.0%),  and 
79.4%.  For  30%  accuracy  at  the  50%  heterogeneity  level  (Table  7-8),  the  CGNN  had  the 
confidence  levels  of  88.5%  for  the  5%  noise  level,  85.6%  for  the  20%  noise  level,  and 
74.4%  for  the  50%  noise  level.  For  the  MLEM,  these  confidence  level  were  91.1%), 
88.8%),  and  76.8%).  At  the  20%  heterogeneity  and  50%  noise  level  (Table  7-7),  the 
confidence  level  of  50%  accuracy  for  the  CGNN  was  93.8%  and  for  the  MLEM  was 
96.8%).  At  the  50%  heterogeneity  and  50%)  noise  level  (Table  7-8),  the  confidence  level 
of  50%  accuracy  for  the  CGNN  was  93.5%  and  for  the  MLEM  was  95.5%. 

When  detailed  density  information  was  not  available  and  only  an  average  density 
was  known,  using  the  2-energy  CGNN  reconstruction  scheme,  the  most  likely  accuracy 
was  50%)  error  with  a  noise  level  as  high  as  20%.  Using  the  2-energy  MLEM 
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Table  7-7.  Comparisons  of  confidence  levels  for  activity  reconstruction  using  perturbed 
kernel  table  and  base  density  kernel  table  under  20%  heterogeneity  condition 


Category  (iv) 

20%  Heterogeneity  with  Local 
Nonuniformity  and  Statistical  Noise 

Error    <50% 

<30% 

<25% 

<20% 

<10% 

C2N05KE2* 
C2N05BE2 

0.998 
0.990 

0.964 
0.911 

0.928 
0.839 

0.857 
0.753 

0.557 
0.468 

C2N20KE2 
C2N20BE2 

0.998 
0.984 

0.951 
0.877 

0.896 
0.810 

0.805 
0.718 

0.495 
0.425 

C2N50KE2 
C2N50BE2 

0.979 
0.938 

0.821 
0.761 

0.743 
0.661 

0.630 
0.543 

0.314 
0.284 

M2N05KE2 
M2N05BE2 

0.998 
0.996 

0.970 
0.939 

0.944 
0.884 

0.866 
0.802 

0.588 
0.499 

M2N20KE2 
M2N20BE2 

0.999 
0.995 

0.965 
0.900 

0.912 
0.846 

0.826 
0.752 

0.524 
0.455 

M2N50KE2 
M2N50BE2 

0.990 
0.968 

0.840 
0.794 

0.757 
0.687 

0.639 
0.571 

0.328 
0.295 

'Notation  for  the  1st  column-  C2:  CGNN  and  20%  heterogeneity; 
M2:  MLEM  and  20%  heterogeneity;  E2:  2-energy; 
K:  perturbed  density  kernel  table;  B:  base  density  kernel  table; 
N05:  5%  noise;  N20:  20%  noise;  N50:  50%  noise. 
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Table  7-8.  Comparisons  of  confidence  levels  for  activity  reconstruction  using  perturbed 
kernel  table  and  base  density  kernel  table  under  50%  heterogeneity  condition 


Category  (v) : 

50%  Heterogeneity  with  Local 
Nonuniformity  and  Statistical  Noise 

Error    <50% 

<30% 

<25% 

<20% 

<10% 

C5N05KE2* 
C5N05BE2 

0.998 
0.984 

0.954 
0.885 

0.918 
0.830 

0.844 
0.741 

0.531 
0.449 

C5N20KE2 
C5N20BE2 

0.998 
0.977 

0.942 
0.856 

0.892 
0.785 

0.789 
0.688 

0.458 
0.406 

C5N50KE2 
C5N50BE2 

0.978 
0.935 

0.818 

0.744 

0.719 
0.669 

0.605 
0.528 

0.315 
0.286 

M5N05KE2 
M5N05BE2 

0.996 
0.993 

0.953 
0.911 

0.916 
0.853 

0.850 
0.761 

0.573 
0.462 

M5N20KE2 
M5N20BE2 

0.997 
0.992 

0.951 
0.888 

0.893 
0.806 

0.814 
0.713 

0.498 
0.408 

M5N50KE2 
M5N50BE2 

0.986 
0.955 

0.828 
0.768 

0.738 
0.675 

0.628 
0.548 

0.319 
0.290 

'Notation  for  the  1st  column-  C5:  CGNN  and  50%  heterogeneity; 
M5:  MLEM  and  50%  heterogeneity;  E2:  2-energy; 
K:  perturbed  density  kernel  table;  B:  base  density  kernel  table; 
N05:  5%  noise;  N20:  20%  noise;  N50:  50%  noise. 
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reconstrruction  scheme,  the  most  likely  accuracy  was  50%  error  with  a  noise  level  as  high 
as  50%. 

Simulations  Using  Eight  Detector  Positions 

All  the  simulations  described  above  were  based  on  four  orthogonal  detector 
positions  90  degrees  apart.  In  this  section,  the  simulation  is  extended  to  eight  detector 
positions  45  degrees  apart. 
Results  of  Activity  Reconstruction  Based  on  the  Perturbed  Kernel  Table 

If  detailed  density  information  is  provided,  total  activity  can  be  reconstructed  by 
using  a  perturbed  kernel  table  which  is  generated  based  on  known  density  distribution. 
Table  7-9  shows  the  simulation  results.  For  both  1 -energy  CGNN  and  1 -energy  MLEM 
reconstructions,  even  with  detector  positions  increased  to  eight,  the  confidence  level  was 
still  poor.  However,  when  increasing  the  energy  to  two,  almost  all  the  predictions  of  both 
the  CGNN  and  the  MLEM  were  within  10%  error  under  the  homogeneous  and  noise-free 
condition. 

For  the  20%  heterogeneity  and  5%  noise  cases,  the  CGNN  predicted  a  20%  error  at  a 
97.2%  confidence  level,  and  the  MLEM  predicted  a  20%)  error  at  a  98.7%)  confidence 
level.  These  results  when  compared  with  the  results  of  four  detector  positions  as  listed  in 
Table  7-4  which  predict  20%  error  at  only  85.7%  confidence  for  the  CGNN  and  only 
86.6%)  confidence  for  the  MLEM,  show  the  great  accuracy  improvement  using  eight 
detector  positions.  The  improvements  were  also  significant  at  the  20%  noise  level.  The 
CGNN  had  a  confidence  level  of  97.9%)  for  25%  accuracy  and  the  MLEM  had  a 
confidence  level  of  95.4%  for  20%  accuracy.  At  the  50%  noise  level,  both  the  CGNN 
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Table  7-9.  Confidence  levels  for  activity  estimation  using  different  reconstruction 
algorithms  and  eight  detector  positions 


Error 

<50% 

<30% 

<25% 

<20% 

<10% 

Homogeneous  and  Noise-Free  Conditions 

HCE1-D8' 
HCE2-D8 

0.676 
1.000 

0.482 
1.000 

0.440 
1.000 

0.391 
1.000 

0.272 
1.000 

HME1-D8 
HME2-D8 

0.747 
1.000 

0.489 
1.000 

0.434 
1.000 

0.372 
1.000 

0.244 
0.996 

20%  Heterogeneitv  and  Noise-Perturbed  Conditions 

C2N05KE2-D8 
C2N20KE2-D8 
C2N50KE2-D8 

1.000 
1.000 
0.994 

0.995 
0.991 
0.894 

0.991 
0.979 
0.814 

0.972 
0.941 
0.715 

0.776 
0.681 
0.404 

M2N05KE2-D8 
M2N20KE2-D8 
M2N50KE2-D8 

1.000 
1.000 
0.990 

0.997 
0.996 
0.883 

0.996 
0.988 
0.802 

0.987 
0.954 
0.685 

0.839 
0.718 
0.381 

50%  Heteroeeneitv  and  Noise-Perturbed  Conditions 

C5N05KE2-D8 
C5N20KE2-D8 
C5N50KE2-D8 

1.000 
0.999 
0.990 

0.994 
0.981 
0.892 

0.983 
0.963 
0.815 

0.955 
0.928 
0.697 

0.757 
0.669 
0.396 

M5N05KE2-D8 
M5N20KE2-D8 
M5N50KE2-D8 

1.000 
0.999 
0.987 

0.996 
0.992 
0.874 

0.994 
0.983 
0.792 

0.978 
0.945 
0.671 

0.822 
0.698 
0.373 

'  Notation  for  the  1st  column-  HC:  CGNN  and  homogeneous  density; 
HM:  MLEM  and  homogeneous  density;  El:  1 -energy;  E2:  2-energy; 
C2:  CGNN  and  20%  heterogeneity;  M2:  MLEM  and  20%  heterogeneity; 
C5:  CGNN  and  50%  heterogeneity;  M5:  MLEM  and  50%  heterogeneity; 
K:  perturbed  density  kernel  table;  D8:  eight  detector  positions; 
N05:  5%  noise;  N20:  20%  noise;  N50:  50%  noise. 
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and  the  MLEM  predicted  30%  error  at  a  confidence  level  lower  than  90%,  but,  they  both 
predicted  50%  error  with  greater  than  99%  confidence.  For  CGNN,  the  most  likely 
accuracy  was  found  to  be  20%  error  with  a  noise  level  up  to  5%.  For  MLEM,  the  most 
likely  accuracy  was  found  to  be  20%  error  with  a  noise  level  up  to  20%.  Figure  7-8 
shows  the  distributions  of  relative  prediction  error  under  20%  heterogeneity  and  different 
noise-perturbed  conditions.  The  maximum  relative  error  for  the  CGNN  at  the  5%,  20% 
and  50%  noise  levels  was  +45%,  +  45%  and  +80%  respectively,  and  for  the  MLEM  was 
+35%,  40%,  and  +80%  respectively.  At  the  20%  and  50%  noise  level,  the  MLEM  tended 
to  underestimate  total  activity,  while  the  CGNN  showed  no  significant  bias. 

In  the  50%)  heterogeneity  condition,  the  overall  performances  of  the  CGNN  and 
the  MLEM  were  only  a  little  worse  than  those  in  the  20%  heterogeneity  condition.  At  the 
5%  noise  level,  the  CGNN  predicted  a  20%  error  at  a  95.5%  confidence  level,  and  the 
MLEM  predicted  a  20%  error  at  a  97.8%  confidence  level.  Compared  to  the  results  of 
the  four  detector  positions  as  listed  in  Table  7-5,  which  predicted  a  20%  error  at  only 
84.4%  confidence  level  for  the  CGNN  and  only  85.0%  confidence  level  for  the  MLEM, 
the  results  of  the  eight  detector  positions  were  greatly  improved.  For  other  noise  levels, 
improvements  were  also  significant.  The  CGNN  predicted  a  25%  error  with  96.3% 
confidence  at  the  20%  noise,  and  50%  error  with  99.0%  confidence  at  the  50%  noise. 
MLEM  predicted  25%  error  with  98.3%  confidence  at  the  20%  noise,  and  50%  error  with 
98.7%  confidence  at  the  50%)  noise.  The  most  likely  accuracy  was  found  to  be  20%  error 
with  a  noise  level  up  to  5%  for  both  the  CGN'N  and  the  MLEM.  Figure  7-9  shows  the 
distributions  of  relative  prediction  error  under  the  50%  heterogeneity  and  different 
noise-perturbed  conditions.  The  maximum  relative  error  for  the  CGNN  at  5%,  20%  and 
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Distributions  of  relative  prediction  error  for  2-energy  CGNN  and  MLEM 
detector  positions)  under  the  50%  heterogeneity  and  different  noise-perturbed 
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50%  noise  level  was  +40%,  +45%  and  +80%  respectively,  and  for  the  MLEM  was  +35%, 
40%,  and  +80%  respectively.  Both  the  CGNN  and  the  MLEM  showed  a  tendency  toward 
underestimating  total  activity  at  the  50%  noise  level. 
Results  of  Activitv  Reconstruction  Rased  on  the  Rase  Densitv  Kernel  Table 

If  detailed  density  information  is  not  provided,  the  reconstruction  of  total  activity 
relies  on  using  the  base  density  kernel  table.  Table  7-10  shows  the  results  of  total  activity 
reconstruction  using  the  base  density  kernel  table  under  heterogeneous  and  5%,  20%,  and 
50%  noise  conditions.  At  the  20%  heterogeneity  condition  with  50%  accuracy,  the 
CGNN  had  confidence  levels  of  98.2%,  97.4%,  and  96.0%,  for  the  three  noise  levels 
respectively;  and  the  MLEM  had  confidence  levels  of  99.0%,  97.9%,  and  95.4% 
respectively.  These  results  were  not  better  than  the  results  of  using  four  detector 
positions  as  shown  in  Table  7-7.  For  50%  accuracy,  the  confidence  levels  shown  in 
Table  7-7  for  the  CGNN  using  base  density  kernel  table  were  99.0%,  98.4%,  and  93.8% 
respectively;  and  for  the  MLEM  were  99.6%,  99.5%,  and  96.8%  respectively. 
Comparing  the  confidence  levels  for  30%  or  25%  accuracy,  using  eight  detector  positions 
were  also  no  better  than  using  four  detector  positions.  These  results  were  the  opposite  of 
the  results  described  in  the  last  two  paragraphs  in  which  the  improvement  of  accuracy  of 
using  eight  detector  positions  was  found.  This  phenomenon  can  be  explained  as  the 
result  of  uncertainty  propagation  in  solving  a  linear  algebra  system.  The  linear  system  of 
2-energy  and  eight  detector  positions  is  16x37,  compared  with  8x37  of  a  2-energy  and 
four  detector  position  system.  Within  a  certain  limit  of  uncertainty,  a  16x37  system  does 
provide  a  more  accurate  prediction  than  a  8x37  system  because  more  equations  are 
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Table  7-10.  Confidence  levels  for  activity  estimation  using  the  base  density  kernel  table 
and  eight  detector  positions 


Error 

<50% 

<30% 

<25% 

<20% 

<10% 

20%  Heterogeneitv  and  Noise-Perturbed  Conditions 

C2N05BE2-D8* 

0.982 

0.905 

0.858 

0.760 

0.490 

C2N20BE2-D8 

0.974 

0.895 

0.841 

0.757 

0.458 

C2N50BE2-D8 

0.960 

0.819 

0.723 

0.619 

0.323 

M2N05BE2-D8 

0.990 

0.888 

0.827 

0.738 

0.471 

M2N20BE2-D8 

0.979 

0.873 

0.810 

0.719 

0.443 

M2N50BE2-D8 

0.954 

0.768 

0.676 

0.568 

0.327 

50%  Heteroseneitv  and  Noise-Perturbed  Conditions 

C5N05BE2-D8 

0.975 

0.882 

0.811 

0.719 

0.444 

C5N20BE2-D8 

0.972 

0.882 

0.801 

0.710 

0.398 

C5N50BE2-D8 

0.954 

0.784 

0.690 

0.575 

0.314 

M5N05BE2-D8 

0.983 

0.854 

0.795 

0.680 

0.399 

M5N20BE2-D8 

0.971 

0.844 

0.775 

0.672 

0.383 

M5N50BE2-D8 

0.946 

0.729 

0.650 

0.546 

0.318 

Notation  for  the  1st  column-  E2:  2-energy; 

C2:  CGNN  and  20%  heterogeneity;  M2:  MLEM  and  20%  heterogeneity; 

C5:  CGNN  and  50%  heterogeneity;  M5:  MLEM  and  50%  heterogeneity; 

B:  base  density  kernel  table;  D8:  eight  detector  positions; 

NQ5:  5%  noise;  N20:  20%  noise;  N50:  50%  noise. 
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provided  to  solve  the  underdetermined  linear  algebra  problem.  However,  if  the 
uncertainty  increases,  as  in  the  situation  of  using  the  base  density  kernel  table  which  adds 
an  extra  noise  component  to  ideal  detector  responses,  the  solution  of  a  larger  linear 
system  may  become  worse  because  more  matrix  manipulations  are  usually  accompanied 
by  more  severe  uncertainty  propagation. 

At  the  50%  heterogeneity  condition  with  50%  accuracy,  the  CGNN  had 
confidence  levels  of  97.5%,  97.2%,  and  95.4%,  for  the  5%,  20%  and  50%  noise  levels 
respectively,  and  the  MLEM  had  confidence  levels  of  98.3%,  97.1%,  and  94.6% 
respectively.  With  30%  accuracy,  the  CGNN  provided  confidence  levels  of  88.2%, 
88.2%,  and  78.4%  respectively,  and  MLEM  provided  confidence  levels  of  85.4%,  84.4%, 
and  72.9%  respectively.  Comparing  these  results  with  the  resuhs  of  using  base  density 
kernel  table  shown  in  Table  7-8,  it  again  shows  that  using  eight  detector  positions  was  not 
better  than  using  four  detector  positions  if  the  base  density  kernel  table  was  used  for 
reconstructing  total  activity  under  heterogeneous  conditions. 

Comparing  the  results  of  the  CGNN  and  the  MLEM  under  the  same  heterogeneity 
and  noise  condition  in  Table  7-10,  it  shows  that  the  MLEM  does  not  always  predict 
higher  confidence  levels  than  the  CGNN.  These  results  were  also  the  opposite  of  the 
results  of  four  detector  positions  as  shown  in  Table  7-7  and  Table  7-8  in  which  the 
superiority  of  using  the  MLEM  over  the  CGNN  is  shown.  This  finding  supports  the  fact 
that  the  MLEM  and  CGNN  reconstruction  algorithms  have  a  limit  of  accurate  estimation 
of  total  activity  depending  on  density  heterogeneity,  noise  level,  and  number  of  detector 
positions.  If  detailed  density  information  is  not  provided  for  a  heterogeneous  barrel  as 
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modeled  in  this  research,  the  use  of  eight  detector  position  does  not  predict  total  activity 
better  than  does  the  use  of  four  detector  positions. 

Simulations  Using  Different  Base  Densities 

The  base  density  (average  density)  used  in  Equation  7-1  is  1.0  g/cml  However, 
the  absorbing  medium  of  most  nuclear  waste  barrels  originated  from  reprocessing  plants 
have  an  average  density  lower  than  1.0  g/cml  This  section  compares  the  total  activity 
estimations  resulted  from  2-energy  CGNN  reconstruction  by  using  different  base 
densities  (0.25,  0.5,  0.75,  1.00  and  1.25  g/cm^)  under  the  50%  density  heterogeneity  and 
20%  noise  condition.  Four  detector  positions  were  used  in  the  simulation. 
Figure  7-10  shows  the  confidence  levels  of  activity  estimations  by  using  the  perturbed 
kernel  table  (the  upper  part  of  the  figure)  and  the  base  density  kernel  table  (the  lower  part 
of  the  figure)  in  the  2-energy  CGNN  reconstruction.  When  base  density  increased  to  1 .25 
g/cm',  confidence  levels  decreased.  However,  confidence  level  of  50%  accuracy  was  still 
above  95%  either  using  the  perturbed  kernel  table  or  using  the  base  density  kernel  table. 
When  base  density  decreased  from  1.00  g/cm\  confidence  levels  increased  and  became 
saturated  at  the  density  lower  than  0.5  g/cm'.  At  0.75  g/cm',  confidence  level  of  30% 
accuracy  was  above  95%  either  using  perturbed  kernel  table  or  using  base  density  kernel 
table. 

Total  Activity  Estimation  for  the  Full  Length  of  the  Waste  Barrel 

All  the  simulations  and  analyses  presented  above  emphasized  on  the  total  activity 
estimation  of  a  vertical  segment  of  the  waste  barrel.  In  the  computer  model,  a  vertical 
segment  is  8  cm  in  height.  For  the  208  ^nuclear  waste  barrel  (88  cm  in  height),  a  total  of 
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Figure  7-10.  Confidence  levels  of  activity  estimations  by  using  perturbed  kernel  table  and 
base  density  kernel  table  in  the  2-energy  CGNN  method  under  different  base  density 
conditions 
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1 1  segments  needs  to  be  analyzed.  This  section  presents  a  statistical  process  to  evaluate 
the  most  likely  accuracy  of  total  activity  estimation  for  the  full  length  of  a  208  ^barrel. 

The  prediction  error  of  total  activity  of  a  vertical  segment  under  a  specific  density 
heterogeneity  and  a  noise-perturbed  condition,  was  characterized  by  a  relative  error 
distribution  function  (see  Figure  7-4  or  7-5).  Statistically,  if  every  segment  has  the  same 
density  heterogeneity  and  noise-perturbed  condition,  the  same  distribution  function  can 
be  used  to  estimate  the  activity  for  each  segment.  Total  activity  for  the  full  length  of  a 
barrel  is  the  sum  of  these  segmented  activities.  To  obtain  the  most  likely  accuracy  of 
total  activity  estimation,  a  Monte  Carlo  sampling  process  was  performed.  First,  activity 
of  each  vertical  segment  was  assumed  to  vary  in  a  range  from  10*  to  10^  Bq  and  was 
assigned  randomly  by  using  the  following  relation: 

^T,  =10^^^''"^     Bq  (7-5) 

where  ST^  is  the  source  activity  for  wth  vertical  segment  and  77^  is  the  random  number 
generated  for  the  mth.  vertical  segment.  Total  activity  {ST)  of  a  barrel  was 


ST  =  Y.ST^ 


(7-6). 


Second,  activity  estimation  error  for  each  vertical  segment  was  sampled  from  a  relative 
error  distribution  function  by  using  a  Monte  Carlo  sampling  technique.  The  estimated 
activity  of  each  vertical  segment  was  obtained  by  adding  the  error  to  the  activity  specified 
from  Equation  7-5.  A  total  activity  was  summed  from  the  estimated  activity  of  each 
segment.  This  estimated  total  activity  was  then  compared  with  ST  to  obtain  a  total 
activity  estimation  error  for  the  full  length  of  a  barrel.  The  above  mentioned  process  was 
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repeated  for  a  total  of  10,000  times  to  define  the  error  bound  and  its  associated 
confidence  level.  Another  type  of  vertical  activity  variation  was  also  studied.  The 
acfivity  of  each  vertical  segment  was  randomly  assigned  in  the  range  of  10'  to  10^  Bq  by: 

5r„  =  10'=-^'-'     Bq  (7-6). 

The  same  Monte  Carlo  sampling  process  described  above  was  used  to  predict  total 
activity  for  the  full  length  of  a  barrel  in  this  activity  range. 

Figure  7-1 1  shows  six  relative  error  distribution  functions  used  in  the  Monte 
Carlo  sampling  process.  These  distribution  functions  were  generated  from  the  2-energy 
CGNN  method  using  the  perturbed  kernel  table  (assuming  a  detailed  density  distribution 
was  provided)  and  the  base  density  kernel  table  (assuming  an  average  density  was 
known)  respectively.  Four  detector  positions  were  used  in  the  simulation.  The  density 
heterogeneity  level  was  50%  and  the  noise  level  was  5%,  20%,  and  50%  respectively. 
The  resulting  confidence  levels  of  total  activity  estimation  for  the  full  length  of  a  208  i 
waste  barrel  are  shown  in  Table  7-12.  For  the  vertical  activity  variation  from  10'  to  10^ 
Bq  and  with  noise  level  as  high  as  50%,  the  most  likely  accuracy  was  20%  error  (98.2% 
confidence  level)  if  a  detailed  density  distribution  was  provided  (case  C5N50KE2),  and 
was  25%.  error  (96.7%  confidence  level)  if  only  an  average  density  was  known  (case 
C5N50BE2).  For  the  vertical  activity  variation  from  10'  to  10^  Bq  and  with  noise  level  as 
high  as  50%,  the  most  likely  accuracy  was  25%  error  (98.1%  confidence  level)  if  a 
detailed  density  distribution  was  provided  (case  C5N50KE2),  and  was  30%  error  (96.5% 
confidence  level)  if  only  an  average  density  was  known  (case  C5N50BE2). 
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Figure  7-11.  Six  relative  error  distribution  functions  used  in  the  total  activity  estimation 
for  the  fiall  length  of  the  waste  barrel 
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Table  7-11.  Confidence  levels  of  activity  estimation  for  the  full  length  of  a  208  ^  waste 
barrel  under  50%  heterogeneity  and  different  noise-perturbed  conditions 


Error   <50%        <30%         <25%         <20%        <10% 


C5N05KE2* 

1.000 

1.000 

0.999 

0.998 

0.921 

C5N20KE2 

1.000 

1.000 

0.999 

0.998 

0.893 

C5N50KE2 

1.000 

0.999 

0.997 

0.982 

0.764 

C5N05BE2 

1.000 

0.999 

0.998 

0.981 

0.733 

C5N20BE2 

1.000 

0.999 

0.992 

0.965 

0.685 

C5N50BE2 

1.000 

0.989 

0.967 

0.914 

0.616 

Vertical  Source 

Strength  Variation  : 

10'~10'Bq 

C5N05KE2 

1.000 

0.998 

0.997 

0.992 

0.852 

C5N20KE2 

1.000 

0.999 

0.998 

0.990 

0.813 

C5N50KE2 

1.000 

0.995 

0.981 

0.940 

0.656 

C5N05BE2 

1.000 

0.996 

0.984 

0.950 

0.684 

C5N20BE2 

1.000 

0.991 

0.973 

0.926 

0.633 

C5N50BE2 

0.998 

0.965 

0.927 

0.859 

0.541 

Notation  for  the  1st  column-  C5:  CGNN  and  50%  heterogeneity;  E2:  2-energy; 
K:  perturbed  density  kernel  table;  B:  base  density  kernel  table; 
N05:  5%  noise;  N20:  20%  noise;  N50:  50%  noise. 


CHAPTER  8 
CONCLUSIONS 

This  research  has  developed  a  new  mathematical  approach  in  finding  an  optimum 
solution  to  a  singular  inverse  radiation  transport  problem  involving  the  nondestructive 
assay  of  radionuclide  inventory  in  a  large  nuclear  waste  container.  The  method 
introduced  is  a  numerical  optimization  scheme  for  determination  of  the  most  likely  value 
of  total  source  activity  based  on  a  pre-specified  maximum  acceptable  error.  First,  a 
computer  model  was  developed,  using  the  CGNN  or  the  MLEM  algorithm  to  reconstruct 
the  source  distribution  of  a  vertical  segment  of  the  waste  barrel  based  on  observed 
external  gamma-ray  detector  responses.  The  computer  model  was  verified  through  36 
experiments  involving  various  source  distribution  patterns  and  various  density 
heterogeneity  conditions  of  the  self-absorbing  medium.  Second,  a  simulation  model  was 
developed  to  consider  a  wide  spectrum  of  measurement  conditions,  including  random 
source  distribution,  density  heterogeneity,  local  nonuniformity  noise,  and  counting 
statistics  noise.  To  simulate  the  density  heterogeneity,  each  computational  grid  was 
randomly  assigned  a  density  perturbed  from  an  average  density  of  1 .0  g/cml  Two 
maximum  density  perturbation  levels,  20%  and  50%,  were  studied.  To  simulate  the  local 
nonuniformity  noise  and  the  counting  statistical  noise,  each  detector  response  was  first 
randomly  perturbed  from  a  noise-free  detector  response  and  then  was  randomly  sampled 
from  a  Gaussian  distribution  fimction.  Noise  perturbation  levels  from  0%  to  50%  were 
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studied.  For  each  density  and  noise  perturbation  condition,  10,000  random  source 
distribution  cases  were  simulated  to  obtain  the  statistical  confidence  level  of  the  specified 
maximum  acceptable  error  for  different  solution  schemes.  The  most  likely  value  of  total 
source  activity  and  the  associated  optimal  solution  scheme  were  thus  determined. 

The  36  experiments  were  performed  by  measuring  the  gamma  signal  of 
distributed  Ir"^  sources  inside  a  208  ^  waste  barrel  with  37  egg-crate  type  compartments 
m  a  vertical  segment.  Measurements  were  made  for  four  orthogonal  detector  positions. 
Total  activity  was  reconstructed  by  the  CGNN  or  the  MLEM  using  detector  responses  of 
two  or  three  energy  peaks.  That  the  computer  model  was  capable  of  accurately 
predicting  total  activity  was  shown  by  the  results  that  all  36  cases  were  within  50%  error 
for  all  the  reconstruction  schemes  used;  32  cases  and  34  cases  were  within  30%  error  for 
2-energy  and  3 -energy  CGNN  reconstruction  respectively;  and  34  cases  and  35  cases 
were  within  30%  error  for  2-energy  and  3-energy  MLEM  reconstruction  respectively. 
The  observed  deviation  of  measured  detector  responses  from  ideal  detector  responses 
suggested  that  local  nonuniformity  noise  and  counting  statistical  noise  need  to  be 
considered  in  the  numerical  simulations. 

Several  findings  are  observed  following  the  numerical  simulations  of  total  activity 
reconstruction  for  a  vertical  segment  of  the  barrel  by  using  four  detector  positions: 

1 .  The  superiority  of  using  multiple  energy  peaks  in  reconstruction  schemes  is 
evident  from  the  simulations  of  the  homogeneous  and  noise-free  condition  which  show 
that  both  the  CGNN  and  MLEM  predict  25%  error  at  100%  confidence  level  using  two  or 


104 

three  energy  peaks,  while  they  predict  25%  error  at  less  than  32.3%  confidence  level 
using  only  one  energy  peak, 

2.  In  noise-free  conditions,  3 -energy  CGNN  reconstruction  enhances  the 
prediction  ability  of  2-energy  CGNN  reconstruction.  However,  in  noise-perturbed 
conditions,  3 -energy  CGNN  reconstruction  deteriorates  the  prediction  ability  of  2-energy 
CGNN  reconstruction.  Compared  with  2-energy  MLEM  reconstruction,  the  prediction 
ability  of  3-energy  MLEM  reconstruction  is  better  at  a  noise  level  lower  than  50%  and 
becomes  worse  when  the  noise  level  is  as  high  as  50%. 

3 .  The  superiority  of  using  the  2-energy  CGNN  or  MLEM  algorithm  over  the 
traditional  SGS  method  is  demonstrated  from  simulations  of  the  homogeneous  and 
noise-free  condition  which  show  that  the  SGS  predicts  25%  error  at  only  27.9%) 
confidence  level  even  when  64  detector  positions  are  used,  while  the  2-energy  CGNN  or 
MLEM  algorithm  predicts  25%  error  at  100%  confidence  level  by  using  only  4  detector 
positions. 

4.  The  confidence  level  of  high  accuracy  becomes  worse  when  the  noise  level  is 
increased.  This  research  incorporates  a  noise  level  as  high  as  50%)  in  the  analyses  to 
represent  an  extreme  condition  of  a  highly  nonuniform  source  distribution  inside  a 
computational  grid.  It  represents  a  high  level  of  detector  response  perturbation  that  may 
result  from  a  point  source  as  small  as  the  Ir''^  piece  used  in  the  experiments.  A  noise 
level  of  5%)  to  10%)  is  a  reasonable  level  for  nuclear  waste  which  usually  has  radioactive 
material  in  a  cluster  and  not  at  a  point.  If  detailed  density  information  is  provided,  total 
activity  can  be  reconstructed  by  using  a  perturbed  kernel  table  which  can  be  generated 
based  on  the  known  density  distribution  of  the  absorbing  medium.  For  the  2-energy 
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CGNN  solution  scheme,  the  most  likely  prediction  error  is  30%  with  a  95.1%  confidence 
level  at  the  20%  heterogeneity  and  20%  noise  condition.  For  the  2-energy  MLEM 
solution  scheme,  the  most  likely  prediction  error  is  30%  with  a  96.5%  confidence  level  at 
the  20%  heterogeneity  and  20%  noise  condition.  At  the  50%  heterogeneity  condition,  the 
most  likely  prediction  error  is  30%  with  a  95.5%  confidence  level  for  the  2-energy 
CGNN  solution  scheme  at  the  10%  noise  level,  and  is  30%  with  a  95.1%  confidence  level 
for  the  2-energy  MLEM  solution  scheme  at  the  20%  noise  level.    Both  the  2-energy 
CGNN  and  MLEM  reconstruction  schemes  predicted  a  50%  error  with  more  than  95% 
confidence  level  at  a  noise  level  as  high  as  50%.  The  high  confidence  level  of  predicting 
error  within  50%  even  at  high  noise  condition  verifies  the  observation  that  36 
experimental  results  are  all  within  50%  error.  If  detailed  density  information  is  not 
provided,  the  reconstruction  of  total  activity  relies  on  using  a  base  density  kernel  table 
which  can  be  generated  based  on  a  known  average  density  of  the  absorbing  medium.  In 
this  situation,  the  most  likely  prediction  error  is  50%  with  more  than  a  95%  confidence 
level  for  both  2-energy  CGNN  and  2-energy  MLEM  solution  schemes  at  the  20%  noise 
level. 

The  computer  model  can  be  extended  to  analyze  a  detection  system  with  more 
than  four  detector  positions.  The  benefit  of  using  eight  detector  position  is  significant  if 
detailed  density  information  is  provided.  For  the  2-energy  CGNN  solution  scheme,  the 
most  likely  prediction  error  is  20%>  with  a  97.2%  confidence  level  at  the  20% 
heterogeneity  and  5%  noise  condition.  For  the  2-energy  MLEM  solution  scheme,  the 
most  likely  prediction  error  is  20%  with  a  98.7%  confidence  level  at  the  20% 
heterogeneity  and  5%  noise  condition.  At  the  50%  heterogeneity  condition,  the  most 
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likely  prediction  error  is  20%  with  a  95.5%  confidence  level  for  the  2-energy  CGNN 
solution  scheme  at  the  5%  noise  level,  and  is  20%  with  a  97.8%  confidence  level  for  the 
2-energy  MLEM  solution  scheme  at  the  5%  noise  level.  At  a  noise  level  as  high  as  20%, 
both  2-energy  CGNN  and  2-energy  MLEM  solution  schemes  have  the  most  likely 
prediction  error  of  25%  with  more  than  95%  confidence.  However,  if  detailed  density 
information  is  not  provided,  using  eight  detector  positions  is  not  better  than  using  four 
detector  positions. 

This  research  also  introduces  a  Monte  Carlo  sampling  process  to  determine  the 
most  likely  accuracy  of  total  activity  estimation  for  the  full  length  of  the  waste  barrel.  In 
this  process,  a  simulated  activity  of  each  vertical  segment  is  randomly  assigned  within  a 
specified  activity  range.  The  activity  prediction  error  for  each  vertical  segment  is 
assumed  to  be  characterized  by  a  relative  error  distribution  tunction  and  is  sampled  from 
this  distribution  function.  An  estimated  activity  of  each  vertical  segment  is  obtained  by 
adding  the  error  to  the  simulated  activity.  A  total  activity  for  the  full  length  of  the  waste 
barrel  is  summed  from  the  estimated  activity  of  each  segment.  For  vertical  activity 
variations  from  10^  to  10'  Bq  and  with  both  density  heterogeneity  and  noise  level  as  high 
as  50%,  the  most  likely  accuracy  is  20%  error  with  a  98.2%)  confidence  level  if  detailed 
density  information  was  provided,  and  is  25%  error  with  a  96.7%  confidence  level  if  only 
an  average  density  is  known.  For  vertical  activity  variations  from  10^  to  10'  Bq  and  with 
both  density  heterogeneity  and  noise  level  as  high  as  50%),  the  most  likely  accuracy  is 
25%)  error  with  a  98.1%  confidence  level  if  a  detailed  density  distribution  is  provided,  and 
is  30%o  error  with  a  96.5%  confidence  level  if  only  an  average  density  is  known. 
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The  simulation  model  developed  in  this  research  can  provide  a  wide  spectrum  of 
experimental  conditions  to  be  simulated.  Several  alternative  solution  schemes  can  be 
tested  and  an  optimal  solution  method  can  be  determined  by  examining  the  statistical 
results  of  activity  prediction  under  different  simulation  conditions.  It  is  especially  useful 
to  quantify  the  error  bound  and  its  associated  confidence  level  of  a  measurement  method 
which  can  not  be  done  by  using  a  limited  number  of  real  experiments. 

Incorporating  nondestructive  measurement  of  multiple  gamma  energy  peaks  in  the 
CGNN  or  MLEM  reconstruction  scheme  is  a  new  concept  developed  by  this  research. 
The  effectiveness  of  this  concept  is  demonstrated.  The  approach  used  in  this  research  can 
provide  the  nuclear  industry  an  accurate  and  confident  method  for  the  nondestructive 
assay  of  radionuclide  inventory  in  a  large  nuclear  waste  container.  The  traditional  SGS 
method,  which  is  proved  to  be  incapable  of  high  accuracy  measurement,  can  be  modified 
to  use  four  detectors  simultaneously  measuring  multiple  gamma  energy  peaks  for  a 
vertical  segment  and  apply  the  reconstruction  schemes  presented  by  this  research.  The 
accuracy  of  prediction  can  definitely  be  increased  to  the  level  described  in  the  above 
paragraphs.  The  recently  developed  TGS  method,  which  usually  can  provide  detailed 
density  information,  can  become  more  accurate  and  less  time  consuming  than  it  is  now,  if 
the  detectors  measuring  emitted  gamma  signals  can  be  fixed  instead  of  rotating  and  the 
concept  of  using  multiple  energy  peaks  in  the  CGNN  or  the  MLEM  algorithm  is  adopted. 


APPENDIX  A 

DERIVATION  OF  CONJUGATE  GRADIENT  SOLUTION  AS  A  LEAST  SQUARED 

SOLUTION  OF  A  QUADRATIC  FUNCTION 

The  solution  of  5  in  the  linear  system  D  =  ^  •  5*  can  be  derived  from  a  least  squared 
solution  of  a  quadratic  function  which  is  defined  as 

e  =  ^-(D  -  KSy  (D  -  KS)  =  ~{D' D  -  D' KS  -  S' K' D  +  S' K' KS) 

To  find  the  minimum  of  e,  we  need  to  differentiate  the  above  equation  with  respect  to  S 
and  let  the  result  to  be  zero.  That  is 

—4  =  (-D'K  -K'D  +  K'KS  +  S'K'K)- 

=  [-(D'-S'K')K-K'(D-KS)]-- 

2 

=  [-(K'(D-KS)y  -K'(D-KS)]■~ 


[-2■K'(D-KS)]■- 


=~K'(D-KS) 
d    e 


d    S 


=  0     means    D=KS.    Therefore,  5"  is  the  least  squared  solution. 
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APPENDIX  B 
DERIVATION  OF  SCALAR  SEARCH  FACTOR  FOR  CONJUGATE  GRADIENT 

ALGORITHM 


To  derive  the  scalar  search  factor  a„  in  the  equation  5'„^,  =  S„+a^d^,  we  need 

to  differentiate  the  quadratic  error  function  e  with  respect  to  a„,  and  then  let  the  result  to 
be  zero.  That  is 

e  =  (D-  KSy  (D-  KS)  =  D' D-  D' KS - S' K' D  +  S' K' KS    , 
and 

^  =  0  =  -D'Kd„-d:K'D  +  d:K'KiS„+aM  +  (S„+aM'K'Kd„ 

=  ~D'Kd„  +S;XKd„  -dlK'D  +  dlK'KS,,  +2a„d:K'Kd„ 

=  -(D-KSjKd„-d:K'(D-KS„)  +  2a„d'„K'Kd„        . 
The  residual  vector  is  defined  as 

r„=K'(D-KSJ       and      r:=iD-KS,yK 
Therefore,  the  above  differential  equation  becomes 
-r:d„-d:r„+2a„d:K'Kd„=0 
d'r 


and  we  have    a.,  - 


d'K'Kd.. 
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APPENDIX  C 
DERIVATION  OF  ORTHOGONALITY  FACTOR  FOR  CONJUGATE  GRADIENT 

ALGORITHM 

The  conjugate  gradient  algorithm  can  be  viewed  as  a  special  case  of  a  general 
Gram-Schmidt  process  for  generating  mutually  conjugate  vectors  dj  from  a  set  of  linearly 
independent  vectors  w /,..., w^  by  the  formulas 

(i,  =  w, 

d^  =  W2  ~t>2\d\ 

d^  -  W3  -b^id^  -b^2^2 

1-1 
The  Gram-Schmidt  orthogonality  requires 

d'.Adj=0  fori^J  (C-2) 

where  ^  is  a  symmetric  and  positive-definite  matrix.  Taking  the  transpose  of  Equation 

C-1  and  then  doing  iimer  product  with  Ad. ,  we  get 

0  =  u'.Adj-bydjAdj    for  i>J 

u'^Ad.      rlK'Kd^ 
^''^=1]M^=1]¥^^        (let  ^  =  i:'i:   and  «,.=,)  (C-3) 

Taking  the  transpose  of  Equation  C-1  and  then  doing  inner  product  with  rj  ,  we  get,  from 

d i^j  ~  0    for  /  <  j    (Property  of  CG  algorithm,  see  Chapter  4  Equation  4-3),  the 
relation 
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Ill 


i-1 


0  =  <n  -  E ^.k^ri  ^      ">;  =  0      for  i  <  j 


k=l 

If  we  select  U-  =  r.  ,  we  have 

rlr.  =  0       for  i  <  j 
From  Chapter  4  Equation  4-5  ,  we  can  derive 


rj,i=^j-ajK'Kdj 


r;r^^,=r;r^-a/K'Kdj 


r'K'Kd,     = 


1   .J. 


or, 


a 


for       j  =  i 
■r/r^     for     j  =  i-l 


/-I 


else 


and  then  from  Equation  C-3  and  C-6,  we  obtain 


*.= 


r'r 


a,_,dl_,K'Kd^_, 


0 


for     j  =  i-\ 
for    j  <  /■  - 1 


(C-4) 


(C-5) 


(C-6) 


(C-7) 


This  relation  shows  that  Z?..  depends  only  on  the  index  i  and  can  be  abbreviated  to  /?,_, . 
Therefore,  the  parameter  bik  in  the  original  Gram-Schmidt  relation  (Equation  C-1)  can  be 
abbreviated  to  have  only  one  index.  We  can  rewrite  Equation  C-1  as 


i-i 


k=\ 

Equation  4-8  of  Chapter  4  is  proved. 


(C-8) 


APPENDIX  D 

PRINTOUT  OF  THE  KERNEL  RESPONSE  TABLE  OF  TWO  ENERGIES  AND 

FOUR  DETECTOR  POSITIONS  FOR  HOMOGENEOUS  BASE  DENSITY  CASE 


Det.  Grid 


Minimum 


1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


7 
8 

1 

2 
2, 
2, 
1 
4 
5, 
7 

12  8 
7, 
5, 
4, 
1 
1 

2, 
2, 
2, 
1, 
1, 
4, 
6 


13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 

25  8 

26  1 

27  8 


28 
29 
30 
31 

32 
33 
34 
35 


.0132E-08 
.8402E-08 
.0580E-08 
,7504E-07 
,2201E-07 
.5124E-07 
,2225E-07 
,7516E-07 
,3683E-07 
,5577E-07 
,2970E-07 
,3869E-07 
,2843E-07 
,5031E-07 
3603E-07 
,2422E-06 
8064E-06 
4599E-06 
9166E-06 
,4757E-06 
8072E-06 
2352E-06 
2533E-06 
1422E-06 
8697E-06 
0759E-05 
7434E-06 
,1815E-06 
,2366E-06 
,3241E-05 
,3619E-05 
,3105E-05 
,3478E-05 
,3077E-05 
.4484E-04 


316  keV 
Maximum 

2.4985E-07 
2.5387E-07 
2.4960E-07 
7.1620E-07 
8.2607E-07 
8.5063E-07 
8.3338E-07 
7.1100E-07 
1.7656E-06 
2.4368E-06 
2.8764E-06 
2.9722E-06 
2.9030E-06 
2.4286E-06 
1.7790E-06 
6.0323E-06 
8.6095E-06 
1.0603E-05 
1.0995E-05 
1.0649E-05 
8.6281E-06 
6.0706E-06 
2.2646E-05 
3.2646E-05 
4.2317E-05 
4.4202E-05 
4.2591E-05 
3.2425E-05 
2.2558E-05 
1.4095E-04 
1.9014E-04 
2.0407E-04 
1.9352E-04 
1.3916E-04 
1.0412E-03 


Average 

1.5049E-07 
1.4952E-07 
1.5084E-07 
3.7479E-07 
4.5759E-07 
4.8979E-07 
4.5624E-07 
3.7196E-07 
9.1157E-07 
1.2272E-06 
1.5590E-06 
1.6792E-06 
1.5429E-06 
1.2326E-06 
9.1518E-07 
2.8962E-06 
4.2222E-06 
5.5543E-06 
6.0387E-06 
5.5196E-06 
4.2123E-06 
2.8825E-06 
1.0353E-05 
1.5153E-05 
2.0703E-05 
2.3584E-05 
2.0839E-05 
1.5163E-05 
1.0486E-05 
6.2318E-05 
8.6448E-05 
1.0106E-04 
8.7565E-05 
6.2370E-05 
3.8865E-04 


588  keV 
Minimum    Maximum     Average 


2.2363E-07 
2.1977E-07 
2.2462E-07 
4.1746E-07 
5.1128E-07 
5.6828E-07 
5.1175E-07 
4.1769E-07 
8.6896E-07 
1.0708E-06 
1.3552E-06 
1.5288E-06 
1.3532E-06 
1.0618E-06 
8.6759E-07 
2.0153E-06 
2.8088E-06 
3.6885E-06 
4.2842E-06 
3.7091E-06 
2.8107E-06 
2.0059E-06 
5.4568E-06 
7.6734E-06 
1.0698E-05 
1.2731E-05 
1.0565E-05 
7.7161E-06 
5.4420E-06 
2.2870E-05 
3.2563E-05 
4.1086E-05 
3.2458E-05 
2.2772E-05 
1.1124E-04 


5.6568E-07 

5.7344E-07 
5.6522E-07 
1.3337E-06 
1.5096E-06 
1.5477E-06 
1.5205E-06 
1.3255E-06 
2.7549E-06 
3.6584E-06 
4.233  lE-06 
4.3566E-06 
4.2673E-06 
3.6484E-06 
2.773  lE-06 
7.5532E-06 
1.0426E-05 
1.2573E-05 
1.2986E-05 
1.2620E-05 
1.0443E-05 
7.5940E-06 
2.2376E-05 
3.1742E-05 
4.0412E-05 
4.2078E-05 
4.0651E-05 
3.1546E-05 
2.2281E-05 
1.0845E-04 
1.4635E-04 
1.5649E-04 
1.4859E-04 
1.0727E-04 
6.4443E-04 


3.7090E-07 
3.6890E-07 
3.7164E-07 
7.7407E-07 
9.1888E-07 
9.7447E-07 
9.1685E-07 
7.693  lE-07 
1.5772E-06 
2.0454E-06 
2.5204E-06 
2.6892E-06 
2.4986E-06 
2.0535E-06 
1.5821E-06 
4.0370E-06 
5.6588E-06 
7.224  lE-06 
7.7871E-06 
7.1856E-06 
5.6490E-06 
4.0191E-06 
1.1403E-05 
1.6285E-05 
2.1707E-05 
2.445  lE-05 
2.1825E-05 
1.6291E-05 
I.1518E-05 
5.2365E-05 
7.2660E-05 
8.4348E-05 
7.3448E-05 
5.2439E-05 
2.6327E-04 
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Det 

.  Grid 

316  keV 

588  keV 

36 

Minimum 

Maximum 

Average 

Minimum 

Maximum 

Average 

T 

1.9622E-04 

1.2040E-03 

5.3356E-04 

1.5056E-04 

7.4381  E-04 

3.5692E-04 

1 

37 

1.4374E-04 

1.0283E-03 

3.8737E-04 

I.1038E-04 

6.3635E-04 

2.6238E-04 

2 

1 

4.3370E-07 

1.7744E-06 

9.1910E-07 

8.6370E-07 

2.7669E-06 

1.5880E-06 

2 

2 

1.2220E-06 

6.0431  E-06 

2.8912E-06 

1.9882E-06 

7.5664E-06 

4.0302E-06 

2 

3 

4.2148E-06 

2.2605E-05 

I.0549E-05 

5.4166E-06 

2.2317E-05 

1.1576E-05 

2 

4 

1.761IE-07 

7.1559E-07 

3.7757E-07 

4.1956E-07 

1.3328E-06 

7.7898E-07 

2 

5 

5.5028E-07 

2.4280E-06 

1.2349E-06 

1.0617E-06 

3.6465E-06 

2.0567E-06 

2 

6 

1.8094E-06 

8.6675E-06 

4.2076E-06 

2.8139E-06 

1.0482E-05 

5.6435E-06 

2 

7 

6.1532E-06 

3.2634E-05 

1.5113E-05 

7.6823E-06 

3.1750E-05 

1.6247E-05 

2 

8 

2.3165E-05 

1.4004E-04 

6.2687E-05 

2.2832E-05 

1.0787E-04 

5.2665E-05 

2 

9 

8.0217E-08 

2.4874E-07 

1.4981E-07 

2.2381E-07 

5.6359E-07 

3.695  lE-07 

2 

10 

2.2145E-07 

8.3464E-07 

4.5910E-07 

5.1021E-07 

1.5225E-06 

9.2176E-07 

2 

11 

7.2430E-07 

2.8943E-06 

1.5495E-06 

1.3465E-06 

4.2561E-06 

2.5070E-06 

2 

12 

2.4635E-06 

1 .0667E-05 

5.5327E-06 

3.6919E-06 

1.2639E-05 

7. 1991  E-06 

2 

13 

8.8087E-06 

4.2766E-05 

2.0653E-05 

1.0635E-05 

4.0800E-05 

2.1666E-05 

2 

14 

3.3324E-05 

1.8999E-04 

8.7454E-05 

3.2321E-05 

1.4606E-04 

7.3370E-05 

2 

15 

1.4247E-04 

1.0369E-03 

3.8750E-04 

1.0953E-04 

6.4179E-04 

2.6252E-04 

2 

16 

7.8309E-08 

2.5427E-07 

1.4771E-07 

2.1955E-07 

5.7417E-07 

3.6534E-07 

2 

17 

2.5146E-07 

8.491  lE-07 

4.8879E-07 

5.6873E-07 

1.5454E-06 

9.7272E-07 

2 

18 

8.3805E-07 

2.9665E-06 

1.6841E-06 

1.5277E-06 

4.3496E-06 

2.6960E-06 

2 

19 

2.9201  E-06 

1.0972E-05 

6.0680E-06 

4.2885E-06 

1.2963E-05 

7.8199E-06 

2 

20 

1.0729E-05 

4.4255E-05 

2.3532E-05 

I.2703E-05 

4.2119E-05 

2.4406E-05 

2 

21 

4.2912E-05 

2.0417E-04 

1.0162E-04 

4.0927E-05 

1.5656E-04 

8.4745E-05 

2 

22 

1.9530E-04 

1.2028E-03 

5.3581E-04 

1.4986E-04 

7.4315E-04 

3.5825E-04 

2 

23 

8.0260E-08 

2.4970E-07 

1.4936E-07 

2.2392E-07 

5.6534E-07 

3.6863E-07 

2 

24 

2.2317E-07 

8.3065E-07 

4.5533E-07 

5.1350E-07 

1.5161E-06 

9.1529E-07 

2 

25 

7.2788E-07 

2.8840E-06 

1.5532E-06 

I.3522E-06 

4.2427E-06 

2.5125E-06 

2 

26 

2.4629E-06 

1.0615E-05 

5.4941  E-06 

3.6917E-06 

I.2588E-05 

7.I593E-06 

2 

27 

8.7715E-06 

4.2647E-05 

2.0916E-05 

1.0591E-05 

4.0699E-05 

2.1899E-05 

2 

28 

3.3325E-05 

1.9239E-04 

8.7138E-05 

3.2326E-05 

1 .4786E-04 

7.3145E-05 

2 

29 

1.4495E-04 

1.0420E-03 

3.8683E-04 

1.1149E-04 

6.4476E-04 

2.6207E-04 

2 

30 

1.7647E-07 

7.1650E-07 

3.7433E-07 

4.203  lE-07 

1.3342E-06 

7.7323E-07 

2 

31 

5.5627E-07 

2.4123E-06 

1.2419E-06 

1.0712E-06 

3.6281E-06 

2.0664E-06 

2 

32 

1.8191E-06 

8.6048E-06 

4.2028E-06 

2.8274E-06 

1.0417E-05 

5.6374E-06 

2 

33 

6.2205E-06 

3.2573E-05 

1.5089E-05 

7.7596E-06 

3.1699E-05 

1.6228E-05 

2 

34 

2.2998E-05 

1.4156E-04 

6.2646E-05 

2.2698E-05 

1 .0888E-04 

5.2593E-05 

2 

35 

4.3100E-07 

1.7753E-06 

9.2684E-07 

8.5934E-07 

2.7680E-06 

1.5987E-06 

2 

36 

1.2318E-06 

6.0661E-06 

2.8770E-06 

2.0021E-06 

7.5896E-06 

4.0123E-06 

2 

37 

4.2556E-06 

2.2605E-05 

1.0593E-05 

5.4642E-06 

2.2326E-05 

1.1620E-05 

3 

1 

1.4309E-04 

1.0405E-03 

3.8495E-04 

1.0997E-04 

6.4399E-04 

2.6104E-04 

3 

2 

1.9561E-04 

1.2019E-03 

5.3079E-04 

1.5006E-04 

7.4270E-04 

3.5522E-04 

3 

3 

1.4266E-04 

1.0133E-03 

3.7753E-04 

1.0967E-04 

6.2772E-04 

2.5663E-04 

3 

4 

2.3036E-05 

1.4021E-04 

6.2780E-05 

2.2721E-05 

1.0799E-04 

5.2713E-05 

3 

5 

3.3465E-05 

1.9081E-04 

8.6544E-05 

3.2438E-05 

1.4662E-04 

7.2709E-05 

3 

6 

4.2841E-05 

2.0423E-04 

1.0177E-04 

4.0863E-05 

1.5659E-04 

8.4858E-05 
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Det. 

Grid 

316  keV 

588  keV 

^ 

Minimum 

Maximum 

Average 

Minimum 

Maximum 

Average 

T 

3.3412E-05 

1.9198E-04 

8.6568E-05 

3.2408E-05 

1.4756E-04 

7.2710E-05 

3 

8 

2.3334E-05 

1.4032E-04 

6.2802E-05 

2.2973E-05 

1.0804E-04 

5.2737E-05 

3 

9 

4.2309E-06 

2.2477E-05 

1.0443E-05 

5.4344E-06 

2.2240E-05 

1.1483E-05 

3 

10 

6.2703E-06 

3.2374E-05 

1.5110E-05 

7.8105E-06 

3.1490E-05 

1.6244E-05 

3 

11 

8.7260E-06 

4.2642E-05 

2.0848E-05 

1.0545E-05 

4.0690E-05 

2.1842E-05 

3 

12 

1.0759E-05 

4.4234E-05 

2.3380E-05 

1.273  lE-05 

4.2101E-05 

2.4263E-05 

3 

13 

8.7407E-06 

4.2612E-05 

2.0903E-05 

1.0561E-05 

4.0668E-05 

2.1894E-05 

3 

14 

6.2456E-06 

3.2663E-05 

1.5369E-05 

7.7774E-06 

3.1757E-05 

1.6479E-05 

3 

15 

4.2143E-06 

2.2417E-05 

1.0456E-05 

5.4161E-06 

2.2167E-05 

1.1497E-05 

3 

16 

1.2225E-06 

6.0382E-06 

2.8910E-06 

1.9893E-06 

7.5620E-06 

4.0287E-06 

3 

17 

1.8088E-06 

8.543  lE-06 

4.2049E-06 

2.8121E-06 

1.0352E-05 

5.6416E-06 

3 

18 

2.4727E-06 

1.0635E-05 

5.5206E-06 

3.7038E-06 

1.2606E-05 

7.1872E-06 

3 

19 

2.9289E-06 

1.0995E-05 

6.0663E-06 

4.2996E-06 

1.2986E-05 

7.8187E-06 

3 

20 

2.4516E-06 

1.0614E-05 

5.5021E-06 

3.6774E-06 

1.2582E-05 

7.1689E-06 

3 

21 

1.8136E-06 

8.6044E-06 

4.2354E-06 

2.8197E-06 

1.0416E-05 

5.6746E-06 

-1 

22 

1.2323E-06 

5.9433E-06 

2.8565E-06 

2.0027E-06 

7.465  lE-06 

3.9903E-06 

3 

23 

4.3451E-07 

1.7778E-06 

9.2862E-07 

8.6508E-07 

2.771 5E-06 

1.6010E-06 

3 

24 

5.4917E-07 

2.4316E-06 

I.2394E-06 

1.0599E-06 

3.6513E-06 

2.0625E-06 

3 

25 

7.2823E-07 

2.8953E-06 

1.5416E-06 

1.3530E-06 

4.2577E-06 

2.4963E-06 

3 

26 

8.3718E-07 

2.9678E-06 

1.6870E-06 

1.5264E-06 

4.3511E-06 

2.6999E-06 

3 

27 

7.2389E-07 

2.8807E-06 

1.5359E-06 

1.3459E-06 

4.2396E-06 

2.4885E-06 

3 

28 

5.5802E-07 

2.4383E-06 

1.2302E-06 

1.0738E-06 

3.6605E-06 

2.0500E-06 

3 

29 

4.3454E-07 

1.7739E-06 

9.1583E-07 

8.6505E-07 

2.7657E-06 

1.5829E-06 

3 

30 

1.7601E-07 

7.1577E-07 

3.7358E-07 

4.1932E-07 

1.3333E-06 

7.7196E-07 

3 

31 

2.2319E-07 

8.3355E-07 

4.5382E-07 

5.1348E-07 

1.5207E-06 

9.1277E-07 

3 

32 

2.5145E-07 

8.4945E-07 

4.9128E-07 

5.6866E-07 

1.5459E-06 

9.7695E-07 

3 

33 

2.2230E-07 

8.3232E-07 

4.5521E-07 

5.1184E-07 

1.5190E-06 

9.1519E-07 

3 

34 

1.7688E-07 

7.1941E-07 

3.7315E-07 

4.2100E-07 

1.3389E-06 

7.7132E-07 

3 

35 

8.0954E-08 

2.4712E-07 

1.5093E-07 

2.2548E-07 

5.605 lE-07 

3.7189E-07 

3 

36 

7.8457E-08 

2.5416E-07 

1.5021E-07 

2.1990E-07 

5.7399E-07 

3.7040E-07 

3 

37 

8.0336E-08 

2.4916E-07 

1.5056E-07 

2.2409E-07 

5.6432E-07 

3.7101E-07 

4 

1 

4.2472E-06 

2.2437E-05 

1.0376E-05 

5.4507E-06 

2.2197E-05 

1.1419E-05 

4 

2 

1.2209E-06 

6.0053E-06 

2.8787E-06 

1.9869E-06 

7.5261E-06 

4.0153E-06 

4 

3 

4.4073E-07 

1.7379E-06 

9.1364E-07 

8.7538E-07 

2.7177E-06 

1.5802E-06 

4 

4 

2.2985E-05 

1.3840E-04 

6.2106E-05 

2.2655E-05 

1.0664E-04 

5.2205E-05 

4 

5 

6.1553E-06 

3.2871E-05 

1.5234E-05 

7.6883E-06 

3.1925E-05 

1.6354E-05 

4 

6 

1.8041E-06 

8.5319E-06 

4.2I25E-06 

2.8066E-06 

1.0340E-05 

5.6481E-06 

4 

7 

5.4946E-07 

2.4300E-06 

1 .2449E-06 

1.0604E-06 

3.6501E-06 

2.0700E-06 

4 

8 

1.7703E-07 

7.1110E-07 

3.7659E-07 

4.2149E-07 

1.3260E-06 

7.7724E-07 

4 

9 

1.4208E-04 

1.0284E-03 

3.8276E-04 

1 .0924E-04 

6.3765E-04 

2.5966E-04 

4 

10 

3.3346E-05 

1.9077E-04 

8.8079E-05 

3.2347E-05 

1 .4657E-04 

7.3813E-05 

4 

11 

8.8561E-06 

4.2349E-05 

2.0787E-05 

1.0681E-05 

4.0433E-05 

2.1790E-05 

4 

12 

2.4739E-06 

I.0550E-05 

5.471 5E-06 

3.7055E-06 

1.2520E-05 

7.1353E-06 

4 

13 

7.3059E-07 

2.8844E-06 

1.5506E-06 

1.3564E-06 

4.2439E-06 

2.5086E-06 

4 

14 

2.2208E-07 

8.3512E-07 

4.5776E-07 

5.1143E-07 

1.5232E-06 

9.193  lE-07 
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Det. 

Grid 

316  keV 

588  keV 

T? 

Minimum 

Maximum 

Average 

Minimum 

Maximum 

Average 

"7 

8.0414E-08 

2.4898E-07 

1 .4922E-07 

2.2426E-07 

5.6405E-07 

3.6835E-07 

4 

16 

1.9679E-04 

1 .2040E-03 

5.341  lE-04 

1.5090E-04 

7.4381E-04 

3.5728E-04 

4 

17 

4.3210E-05 

2.0387E-04 

1.0215E-04 

4.1180E-05 

1.5636E-04 

8.5138E-05 

4 

18 

1.0780E-05 

4.4277E-05 

2.3568E-05 

1.2754E-05 

4.2137E-05 

2.4430E-05 

4 

19 

2.9203E-06 

1.0978E-05 

6.0576E-06 

4.2892E-06 

1.2969E-05 

7.8087E-06 

4 

20 

8.4088E-07 

2.9697E-06 

1.6741E-06 

1.5319E-06 

4.3534E-06 

2.6821E-06 

4 

21 

2.5115E-07 

8.5022E-07 

4.908  lE-07 

5.6813E-07 

1.5470E-06 

9.7612E-07 

4 

22 

7.8587E-08 

2.5444E-07 

1.4870E-07 

2.201 8E-07 

5.7451E-07 

3.6733E-07 

4 

23 

1.4627E-04 

1.0471E-03 

3.8804E-04 

1.1201E-04 

6.4809E-04 

2.6299E-04 

4 

24 

3.3640E-05 

1.9192E-04 

8.6210E-05 

3.2616E-05 

1.4752E-04 

7.2465E-05 

4 

25 

8.7433E-06 

4.1875E-05 

2.0910E-05 

1.0563E-05 

4.0069E-05 

2.1901E-05 

4 

26 

2.4567E-06 

1.0631E-05 

5.4898E-06 

3.6838E-06 

1.2603E-05 

7.1550E-06 

4 

27 

7.2664E-07 

2.8776E-06 

1.5417E-06 

1.3504E-06 

4.2345E-06 

2.4970E-06 

4 

28 

2.2225E-07 

8.3428E-07 

4.5401E-07 

5.1172E-07 

1.5219E-06 

9.1333E-07 

4 

29 

8.0500E-08 

2.4986E-07 

1.5038E-07 

2.2446E-07 

5.6567E-07 

3.7068E-07 

4 

30 

2.3112E-05 

1.3917E-04 

6.2940E-05 

2.2764E-05 

1.0698E-04 

5.2834E-05 

4 

31 

6.1412E-06 

3.2929E-05 

1.5091E-05 

7.6717E-06 

3.1991E-05 

1.6223E-05 

4 

32 

1.8029E-06 

8.5408E-06 

4.1832E-06 

2.8053E-06 

1.0350E-05 

5.6180E-06 

4 

33 

5.5184E-07 

2.4387E-06 

1.2437E-06 

1.0642E-06 

3.661  lE-06 

2.0681E-06 

4 

34 

1.7730E-07 

7.2015E-07 

3.7514E-07 

4.2173E-07 

1.3401E-06 

7.7481E-07 

4 

35 

4.2227E-06 

2.2310E-05 

1.0379E-05 

5.426  lE-06 

2.2078E-05 

1.1418E-05 

4 

36 

1.2276E-06 

5.9304E-06 

2.8914E-06 

1.9965E-06 

7.444  lE-06 

4.0321E-06 

4 

37 

4.3292E-07 

1.7640E-06 

9.2020E-07 

8.6242E-07 

2.7537E-06 

1.5892E-06 
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